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CHAPTER  I 


INTRODUCTION 

The  current- induced  transition  to  the  normal  state  of  one-dimen¬ 
sional  superconducting  filaments  (ODSF)  is  one  of  the  most  interesting 
problems  in  nonequilibriun  superconductivity.  The  term  one-dimensional 
is  used  for  samples  with  transverse  dimensions  small  with  respect 
to  the  Ginzburg-Landau  (GL)  coherence  length  ^  (T)  and  to  the  London 
penetration  length  X(T),  and  with  length  large  compared  with  those 
parameters.  In  this  situation,  any  spatial  dependence  of  the  supercon¬ 
ducting  properties  of  the  sample  is  assumed  to  occur  only  along  th° 
filament. 

Experimentally,  the  resistive  transition  of  a  superconducting  fila¬ 
ment  can  be  studied  when  a  constant  current  is  applied  to  the  sample 
and  the  voltage  across  the  sample  is  measured  as  a  function  of  the  tem¬ 
perature  (V  versus  T  curves,  also  called  R  vs.  T  curves).  Conversely, 
the  temperature  can  be  kept  fixed  and  the  voltage  across  the  sample 
measured  as  a  function  of  the  current  applied  (V  versus  I  curves,  com¬ 
monly  called  I-V  curves) . 

For  temperatures  close,  but  not  too  close,  to  the  "mean-field" 
critical  temperature  T0  of  the  sample,  that  is,  for  0.1  £ 1-T/T„^ 0.001 , 
the  onset  of  resistance  is  characterized  by  the  appearance  of  reproduc¬ 
ible  regular  voltage  steps  occurring  at  specific  values  of  temperature 
and  current.  This  structure  was  first  observed  in  1968  by  Webb  and 

Warburton1  in  the  R  vs.  T  curve  of  whisker  crystals  of  tin  with  cross- 

o 

sectional  areas  smaller  than  1  (^Mm)  . 

2 

A  few  years  later,  Meyer  and  Minnigerode  published  their  expen- 


1 
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mental  results  on  the  resistive  transition  of  tin  whiskers.  They  ob¬ 
served  voltage  steps  in  the  V  vs.  T  curves  for  current  intensities 
larger  than  35  A  and  found  that  the  height  of  the  voltage  steps  in¬ 
creased  linearly  with  the  current  applied.  They  also  measured  the  I-V 
curves  at  fixed  temperatures  and  observed  voltage  steps  followed  by  a 
plateau  of  constant  differential  resistance.  Figure  1.1  3hows  a  set  of 
I-V  curves  taken  from  reference  2.  These  are  the  first  published  set  of 
I-V  curves  showing  the  step-like  structure  and  for  that  reason  were 
chosen  to  be  the  first  figure  in  this  thesis. 

Similar  I-V  curves  have  been  measured  in  long  microbridges  of  tin^, 
aluminum4,  indiun^  and  niobium^,  and  also  in  indiun  whiskers^. 

Two  very  important  features  of  those  I-V  curves  are,  firstly,  the 
changes  in  differential  resistance  at  all  the  steps  in  a  given  curve  have 
the  same  magnitude,  and  secondly,  the  presence  of  an  excess  current, 
meaning  that  the  extrapolation  to  zero  voltage  of  the  linear  part  of  the 
steps  intercepts  the  current  axis  at  a  finite  positive  value,  Iexc- 

Skocpol,  Beasley  and  Tinkham^  developed  a  model,  known  as  the  SBT 
model,  to  explain  those  features  in  the  I-V  curves  of  superconducting 
filaments.  According  to  SBT,  each  voltage  step  corresponds  to  the  appear¬ 
ance  of  a  spatially  localized  "phase-slip  center"  (PSC)  which,  in  simple 
terms,  can  be  described  as  a  dissipative  unit  that  is  able  to  carry  a  fi¬ 
nite  supercurrent.  Each  unit  consists  of  a  region  of  length  b;  £ (T) 

where  the  G l  order  parameter  undergoes  a  relaxation  oscillation  at  a 
Josephson  frequency  determined  by  the  voltage  across  the  PSC.  This  volt¬ 
age  is  determined  by  the  relaxation  of  the  nonequilibriun  potential 


Figure  I. 1 .(Xir^ent-voltage  characteristics  of  a  tin  whisker  at  various 
temperatures  obtained  by  Meyer  and  Minnlgerode  (reference  2) . 
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which  occurs  in  a  longer  length  scale  corresponding  to  the  diffusion 
length  of  the  quasiparticle  charge,  namely,  -A- ga¬ 
in  order  to  explain  the  significance  of  the  SET  model,  it  is  neces¬ 
sary  to  provide  a  historical  review.  For  that  reason,  the  next  chapter 
is  organized  chronologically.  Starting  with  the  origins  of  the 
concept  of  phase  slippage,  I  will  discuss  briefly  the  models  of  PSCs 
proposed  by  different  authors  before  and  after  SBT.  A  general  view  of 
the  most  relevant  experiments  on  PSCs  will  be  also  presented. 

As  it  will  be  shown  in  the  next  chapter,  all  the  models  agree  in 
the  explanation  of  the  voltage  steps  in  the  I-V  curves  of  ODSF  in  terms 
of  the  appearance  of  PSCs.  However,  the  interpretation  of  the  quasipar¬ 
ticle  diffusion  length  has  been  a  subject  of  controversy  over 

the  last  few  years.  For  that  reason  we  have  performed  a  series  of 

O 

experiments0  in  order  to  measure  the  decay  with  distance  from  the 
center  of  the  PSC  of  the  quasiparticle  current  and  the  nonequilibriim 
potential  and  we  have  compared  our  results  with  the  existing  theories 
and  with  the  experimental  results  obtained  by  other  authors. 

In  one  experiment  we  have  measured  the  effect  that  the  quasiparti¬ 
cle  current  generated  in  one  PSC  has  on  the  appearance  of  other  PSCs  in 
the  same  microbridge.  The  geometry  of  these  samples  was  such  that  the 
nucleation  of  PSCs  in  different  locations  was  controlled  by  applying 
current  through  side  probes  to  differents  segments  of  the  microbridge 
(I  will  refer  to  these  samples  as  samples  of  type  S) .  In  that  way  the 
interaction  effects  could  be  measured  as  a  function  of  the  distance 
between  the  two  interacting  PSCs.  We  have  separated  heating  effects  from 
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diffusive  current  effects  in  the  same  fashion  as  Jillie^  did  is  his 
work  on  interacting  short  tnicrobridges.  We  have  observed  that  the  quasi¬ 
particle  diffusion  current  decays  exponentially  with  distance  and  that 
for  a  temperature  T=  0 .99TC ,  the  magnitude  of  the  characteristic  decay 
length  is  -Aq#=  (6  t  1)yWm,  corresponding  to  a  relaxation  time  of 
“Sq#2  (7  t  2)x10-10  sec  for  our  tin  microbridges.  The  temperature 
rise  of  the  microbridge  caused  by  the  PSC  corresponds  to  a  heat  transfer 
coefficient  =  1.5  watt /K  cm  for  our  tin  sanples  on  glass  substrates 
and  in  direct  contact  with  the  liquid  He  bath.  No  voltage  locking  was 
observed,  presumably  due  to  the  relatively  long  distance  between  PSCs 
compared  to  the  decay  length  of  the  a.c.  nonequilibrium  currents  -A 

3  .C  . 

(Aa_c2:0.6  -Aq,  for  tin  at  T=  0.99TC)10. 

In  a  different  type  of  experiment  we  have  measured  directly  the 
exponential  decay  of  the  quasiparticle  potential.  We  have  fabricated 
tin  micro'oridges  with  normal  tunnel  junction  probes,  similar  to  the 
samples  used  in  the  Dolan  and  Jackel11  experiment  (I  will  refer  to 
these  samples  as  samples  of  type  T) .  We  have  proved  that  the  decay  of 
the  quasiparticle  potential  is,  in  fact,  exponential  and  that  the  char¬ 
acteristic  decay  length  Aq#  diverges  at  the  critical  temperature  as 

A-1  /4 

gfoC  (1  -  T/Tc)  which  corresponds  to  a  relaxation  time 
'Cq,(T)=  (0.7±  0.2)x10“10sec(1  -  T/Tc)_1/2,  which  agrees  very 
well  with  the  result  of  our  previous  experiment. 

We  have  also  injected  normal  current  in  the  microbridge  through 
the  normal  probes  in  samples  of  type  T.  The  induced  nonequilibriun  po¬ 
tential  was  measured  at  other  points  of  the  bridge  and  we  observed  the 
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decay  of  this  potential  with  the  distance  from  the  injection  site.  The 
effect  of  that  injected  non equilibrium  current  on  the  appearance  of  the 
first  PSC  in  the  microbridge  was  also  measured  and  we  observed  the 
decay  of  this  effect  with  the  distance  between  the  injection  site  and 
the  PSC.  The  magnitude  of  the  decay  length  obtained  in  these  experiments 

*  Q 

also  agrees  with  the  results  of  the  other  experiments' . 

We  have  not  found  any  conclusive  evidence  of  the  existence  of  the 

12 

"matching  region" ,  proposed  by  Baratoff  ,  where  the  decay  of  the 
nonequilibrium  potential  is  non-exponential. 

According  to  the  theory  J  very  close  to  Tc,  Gq*  can 

be  written  as  "Sq*=  0.73  'C  E  A  CO)/  A  (T).  Using  this  relation,  we 
have  found  that  the  inelastic  phonon  scattering  time  in  our  tin 

samples  at  Tfi  is  7>E=  (1.6  t  0.4)x10“ 10sec. 

The  last  part  of  this  work  consists  of  the  detection  of  the 
oscillations  of  the  supercurrent  in  the  core  of  the  PSC.  We  have  fabri¬ 
cated  tin  microbridges  of  various  lengths,  from  short  microbridges 
(L  2r  1  pm)  to  long  microbridges  (L~  4 0pm)  and  the  effect  of  an  external 
microwave  field  on  the  I-V  curves  of  those  microbridges  was  investigated. 
We  have  observed  well  defined  Josephson  steps  in  microbridges  up  to  lOy^m 
long. 


CHAPTER  II 

THEORIES  AND  EXPERIMENTS  ON  PHASE-SLIP  CENTERS. 

II. 1  DEFINITION  OF  PHASE-SLIPPAGE 

The  concept  of  a  phase-slip  event  was  developed  by  Langer  and 
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Ambegaokar  (LA)  and  it  is  based  on  an  idea  proposed  by  Little  '. 
According  to  LA,  if  V12  is  the  voltage  difference  between  points  x1 
and  Xp  at  each  er>d  °?  the  filanent  and 

^(x)  =  JyCx)|  ei(P(x)  (II. 1) 

is  the  GL  order  parameter,  then  the  phase  difference  between  those  two 
points  evolves  in  time  following  the  Josephson  relation 

—  (Cp(x2).(p(Xl))=  ^12  (II.  2) 

dt  -h 

According  to  the  GL  theory,  the  magnitude  of  the  supercurrent  in  the 
filament  is 

<n.3> 

In  a  one  dimensional  filament,  the  magnetic  field  can  be  neglected  and 
we  assume  \=0.  Then 

V»>  *  iSM2  a?  (ii. «) 

From  equations  (II. 2)  and  (II. 4)  it  follows  that  in  the  presence  of 
a  voltage  difference,  the  supercurrent  along  the  filament  will  increase 
continually  without  limit.  Therefore,  in  order  to  achieve  a  steady  state 
situation,  thermal  fluctuations  in  the  interior  of  the  superconductor 
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must  occur  in  order  to  reduce  the  phase  difference  at  the  same  rate  as 

the  voltage  increases  it.  According  to  LA,  when  the  supercurrent  J 

reaches  a  critical  value  Jc,  the  phase  coherence  is  broken  at  a  certain 

point  in  the  filament,  and  the  GL  order  parameter  vanishes  briefly, 

allowing  the  phase  to  slip  by  2 IT  before  the  coherence  is  reestablished. 

LA  found  that  fluctuations  of  the  order  parameter  that  allow  phase 

slippage  are  most  likely  to  occur  over  a  length  £  (T) ,  the  GL  coherence 

length.  A  greater  or  lesser  length  is  energetically  less  favorable.  The 

1  ft 

LA  theory,  with  the  correction  introduced  by  McCunber  and  Halperin  , 
describes  very  well  the  onset  of  resistance  in  one-dimensional  supercon¬ 
ducting  filaments  (ODSFs)  in  the  limit  of  low  current  density  and  tem¬ 
peratures  very  close  to  T  ,  i.  e.  1-T/T  4  o.OOI. 

c  c 
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II. 2  PRE-SBT  MODELS  OF  PHASE-SLIP  CENTERS 

The  first  model  intended  to  explain  the  voltage  steps  observed  in 

iq  pn 

the  resistive  transition  of  ODSFs  was  proposed  by  Fink  ’  in  1973  and 
for  the  first  time,  the  expression  "phase-slip  center"  was  used. 

According  to  Fink,  when  the  total  current  along  the  filament 
exceeds  the  critical  current,  the  order  parameter  becomes  zero  at  some 
points  in  the  filament,  permitting  the  phase  of  the  order  parameter  to 
change  in  time  at  those  points.  In  this  model,  the  modulus  of  the  order 
parameter  is  time-independent  but  spatially  periodic.  The  phase  is  a 
time-dependent  quantity  and  its  evolution  is  described  in  terms  of  the 
time  dependent  Ginzburg-Landau  (TDGL)  theory. 

The  voltage  across  each  PSC  is  an  integer  fraction  of  the  total 
voltage  across  the  filament.  Fink  calculated  the  I-V  curves  and 

P 

compared  them  to  the  experimental  results  of  Meyer  et  al.  .  The 
qualitative  agreement  is  only  fair,  the  main  discrepancy  being  related 
to  the  existence  of  an  excess  current.  The  I-V  curves  calculated  by 
Fink  extrapolate  to  the  origin  and  he  suggested  that  the  offset  in  the 
experimental  I-V  curves  was  caused  by  a  recorder  drift. 

Skocpol  et  al.^  also  tried,  unsuccessfully,  to  fit  their  data  to 
Fink's  model. 
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II. 3  SBT  MODEL  OF  A  PHASE-SLIP  CENTER 

•j 

According  to  Skocpol, Beasley  and  Tinkham  (SBT)  ,  when  the  current 
along  a  ODSF  exceeds  the  critical  current  of  the  weakest  point,  the  fil¬ 
ament  jumps  into  a  dissipative  state  with  periodic  phase-slip  events 
occurring  in  that  weakest  point  at  a  Josephson  frequency  set  by  the  total 
voltage  between  the  ends  of  the  filament. 

SBT  considered  a  two-fluid  picture  to  describe  the  current  through 
the  phase-slip  center.  One  part  of  the  current  is  carried  by  quasipar¬ 
ticles,  the  dynamics  of  which  is  determined  by  a  normal  electrochemical 
potential  J\  \  the  other  part  of  the  current  is  carried  by  the  condensate 
with  which  a  superconducting  electrochemical  potential  is  associated. 

Using  the  results  of  the  GL  theory,  SBT  wrote  the  supercurrent 
density  as 


and  from  the  Josephson  relation 

*  =  -2  M 

dt  ‘  s 


(II .4) 


(II. 5) 


Combining  those  two  equations,  one  obtains 


(II. 6) 


An  important  SBT  assumption  is  that  the  normal  current  density  is 
given  by 


where 


c 


—  l£n  (II. 7) 

pe  dx 

is  approximately  the  normal  resistivity  of  the  filament. 
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As  we  mentioned  before,  in  order  for  a  phase-slip  event  to  occur, 
the  order  parameter  (x)|  must  vanish.  The  difference  with  Fink's 
model  is  that,  according  to  SBT,  |  ^l(x)|  vanishes  only  momentarily,  driv 
ing  to  zero,  with  the  total  current  being  carried  by  normal  electrons 
before  ( x) ]  regains  its  maximum  value  to  collapse  again  in  the  next 
phase-slip  event.  This  process  is  just  a  relaxation  oscillation  of  the 

pi 

supercurrent  like  in  the  model  of  Notarys  and  Mercereau  of  the  dy¬ 
namic  behaviour  of  non-tunneling  type  weak  links. 

SBT  calculated  the  time  average  of  the  supercurrent  at  the  center 
of  the  PSC  by  using  current-carrying  solutions  of  the  GL  equations,  so 
that  equation  (II. 4)  can  be  written  in  the  form 

Js=  (1  -  HI- 8) 


where  q=d(p/dx.  When  the  average  is  taking  between  q=0  and  q=1/£  ,the 
result  is 


=0.65  Jc 


where  JQ  is  the  critical  current  density  Jq=2  Jq/3|3. 

According  to  the  SBT  model,  the  oscillations  of  the  order  paraneter 
occur  in  a  length  of  the  order  of  £  (T)  (this  region  is  called  the  core 
of  the  PSC)  as  proposed  by  Langer  and  Ambegaokar1^.  However, the  conver¬ 
sion  of  normal  current  into  supercurrent  occurs  in  a  length  scale  which 
is  determined  by  the  quasiparticle  diffusion  length  A q#, usually  much 

larger  than  £  (T).  This  region  is  called  the  wings  of  the  PSC  in  the 

22 

most  recent  literature  .  The  quasiparticle  diffusion  length  can  be 
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written  as 

-A.q#  =  (II. 9) 

where  Vp  is  the  Fermi  velocity,  l  is  the  mean  free  path  and  is 

the  relaxation  time  for  the  processes  which  lead  to  the  relaxation  of 
the  difference  ^ .  The  interpretation  of  the  time  '"Cq,  will 
be  discussed  in  section  II. 5. 

S8T  also  calculated  the  spatial  dependence  of  the  time  average  of 
the  pair  potential  and  the  normal  potential  They  found  that 

is  constant  at  each  side  of  the  core  and  it  changes  discontinu- 
osly  at  the  center  of  the  PSC.  On  the  other  hand,  y*s  changes 

more  gradually 

) ) s inh/?r^o\  /cosh/^,l,rx,a\ 

\\J/  CII- 10) 

xQ  is  a  point  where  like  at  the  ends  of  the  filament, 

where  it  makes  contact  with  the  banks,  and  is  the  position  of  the 
center  of  the  PSC  where  J^x^):  0-65  J  .  For  a  very  long  filament 
(L  >">  Ag#),with  bias  current  going  from  left  to  right,  equation  (11.10) 
reduces  to 

/V/V  (11.11) 

LV^q./J 

for  points  to  the  left  of  the  core  x  <  x1 ,  and 

/V/V  e^Q»(J-^(x-,))exp  (11.12) 

for  points  to  the  right  of  the  core  x  >  x.j. 


/V/V  -ecV(J-Vxi 
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Fran  (1.7),  the  time  average  of  the  normal  current  is 

(11.13) 


Jn(x)  =  (J-J3(Xl))exp^-lirf.^ 


Figure  II. 1  shows  the  time  average  of  the  electrochemical  potentials 
and  currents  as  a  function  of  the  coordinate  along  the  filament  (from 
reference  3) •  The  total  voltage  drop  across  the  filament  is: 


-w 

wd 


(11.14) 


and  the  differential  resistance  associated  with  the  PSC  is 

(11.15) 


dl 


wd 


where  wi  is  the  cross-sectional  area  of  the  filament.  Equation  (11.15) 
is  one  of  the  most  important  results  of  the  SBT  model.  It  equates 
the  differential  resistance  of  a  PSC  with  the  resistance  of  a  normal 
segment  of  length  2Aq». 

Insofar  as  these  results  of  the  simple  model  are  correct,  the 
magnitude  of  Aq#  and  the  value  of  Ig(x.j)  can  be  inferred  directly 
from  the  experimental  I-V  curves. 

Now,  I  will  discuss  briefly  some  of  the  models  of  PSCs  developed  by 


other  authors. 


I 
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Figure  II. 1. Time  averages  of  (a)  electrochemical  potentials  and 

(b)  supercurrent  and  normal  current  in  a  long  microbnage 
with  a  single  phase-slip  center  according  to  SBT  (reference  _) 


II. 4  POST-SBT  MODELS  OF  PHASE-SLIP  CENTERS 


In  order  to  deal  with  the  time-dependent  character  of  the  PSC, 
Kadin,  Smith  and  Skocpol  (KSS)^  developed  a  clever  model  in  which  a  PSC 
can  be  visualized  as  a  Josephson  oscillator,  of  zero  length  and 
negligible  resistance  connected  to  a  transmission  line.  If  the  total 
current  through  the  PSC  is  constant,  the  normal  current  injected  in  the 
transmission  line  will  have  an  a.c.  character  because  of  the  a.c. 
supercurrent  generated  in  the  Josephson  element.  This  a.c.  normal 
current  decays  more  rapidly  than  its  d.c.  component.  Figure  II. 2  shows 
the  spatial  dependence  of  the  normal  current  at  a  certain  time.  The 
time  average  of  the  current  and  potentials  reduce  to  the  results  of  the 
SBT  model. 

KSS10  also  found  that  the  differential  resistance  of  the  PSC  is 
given  by  equation  (11.15)  and  that  the  extrapolated  excess  current  Iexc 
is 

Jexc  =Ts(x1)  8  <f>  h  (H.16) 


where  <f>  is  the  average  of  the  current-phase  relation  of  the  Josephson 
element 


1  f2TT 

<f>=  JL  f(A0)d(^S) 
27T  L 


Is=  Icf(69(t)) 


(11.17) 


A0  is  the  phase  difference  across  the  core.  Equations (II.  16)  and  (iI.Vj) 
are  valid  if  the  voltage  across  the  Josephson  element  is  constant  so 
that  the  phase  difference  AS  increases  linearly  in  time.  This 
assunption  is  correct  at  the  high  voltage  limit  of  the  I-V  curve  where 
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the  oscillations  of  the  normal  current  are  small  compared  to  the  total 
current.  For  a  linear  current-phase  relation,  equations (II. 16)  and  (II. 17) 

«ive  WV2- 

A  different  approach  has  been  taken  by  other  authors. 
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In  a  series  of  papers  '  '  ,  Galaiko  has  developed  a  model  based 

on  the  microscopic  theory  for  the  resistive  state  of  ODSFs  in  the  range 
of  current  <  j  <  jc2  and  very  close  to  T  .  j  is  the  critical 
current  determined  by  the  dependence  of  the  number  of  superconducting 
electrons  on  the  velocity  of  the  condensate  (pair  breaking  current)  and 
jc2  is  the  current  below  which  superconductivity  pairing  develops  in  the 
normal  state.  In  Galaiko' s  model,  the  resistive  state  of  a 
superconducting  channel  consists  of  a  dynamic  =  rrqcture  cf  normal  and 
superconducting  regions  alternating  along  try*  cnannel,  like  In  Fink's 
model  7'  .In  the  presence  of  an  electric  field,  this  state  is 
characterized  by  discontinuities  of  the  phase  of  the  complex  order 
parameter  in  the  regions  where  &  =0.  Galaiko  calls  those  regions 
"centers  of  singularity".  Galaiko  solved  the  kinetic  equations,  derived 
in  a  previous  paper  0  based  on  the  BCS  theory,  for  a  dirty 
superconductor  of  infinite  length  and  for  temperatures  very  close  to  Tc 
where  the  Joule  heating  can  be  neglected.  The  spatial  period  of  the 
structure  of  normal  and  superconducting  regions  in  the  filament  was 
found  by  maximizing  L  for  a  given  current. 

As  j  increases  approaching  jc2  ,  the  spatial  period  increases 
because  of  the  size  of  the  normal  regions  increasing.  As  j  decreases, 
the  spatial  period  decreases  until  the  superconducting  regions  reach  a 
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size  £(T)  and  then  increases  again  because  of  the  increase  of  the  size 
of  the  superconducting  regions. 

As  j  gets  closer  to  j  ,  the  resistive  state  is  determined  by  the 
penetration  of  the  electric  field  g  from  the  normal  regions  into  the 
superconducting  regions.  Galaiko  calculated  the  shape  of  the  I-V  curve 
with  j  in  the  vicinity  of  and  jc2  and  an  excess  current  was  found  in 
both  limits.  For  j~jQ  the  excess  current  is  associated  with  the 
addition  of  a  superconducting  current  while  at  j £ jc2  it  is  associated 
with  the  increase  of  the  normal  conductivity.  Galaiko' s  model  is  valid 
only  for  infinite  samples;  however,  he  suggests  that  for  samples  of 
finite  length,  "centers  of  singularity"  should  appear  and  that  the 
voltage  steps  observed  experimentally  in  the  1-V  curves  of  ODSFs  are 
related  to  the  penetration  of  new  centers  into  the  sample,  but  no 
numerical  results  for  the  magnitude  of  those  jumps  are  derived  in  his 
papers. 

Kramer  and  Baratoff^  have  solved  numerically  the  time-dependent 

Ginzburg-Landau  (TDGL)  equations  for  homogeneous  ODSFs.  They  found  that 

for  a  limited  range  of  current  jm^n  <  j  <  jc,  localized  phase-slip 

oscillations  occur  in  which  |ip|  goes  to  zero  at  one  point  and  the  phase 

changes  by  2T1  .In  Kraner's  nomenclature,  the  experimental  critical 

current  at  which  a  voltage  first  appears  is  jmax  which  corresponds  to 

the  GL  value  for  J  used  in  the  SBT  model,  j  .  is  the  current  below 
c  -min 

which  the  state  is  pure  superconducting,  j  is  the  critical  current 
above  which  the  superconducting  state  becomes  globally  unstable  or  below 
which  the  normal  state  is  globally  unstable.  As  j  decreases  toward 
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the  period  of  the  oscillation  goes  to  infinity.  Above  j  the  order 

parameter  performs  oscillations  with  decreasing  amplitude  between 

phase-slip  events.  Eventually,  an  expanding  normal  domain 

develops,  i  =  0.385,  j  =  0.335  and  1.  =0.326  in  units  of  j  =  hcV8enX  !• 
iQdx  c  min  o 

28 

Later,  Kramer  and  Watts-Tobin  found  that  phase-slip  oscillations 

are  possible  for  a  larger  range  of  currents.  They  solved  generalized 

29 

TDGL  equations  derived  from  microscopic  theory  and  taking  into  account 
the  effect  of  a  nonequilibrium  quasiparticle  distribution. 

For  y  — »  0,  where  ^  -  2'^Aq/Ti,  their  equations  become  the  simple 
TDGL  equations  derived  in  reference  27.  is  the  inelastic  scattering 

time  and  A0  is  the  temperature-dependent  value  of  the  equilibrium 
order  parameter.  Kraner  and  Watts-Tobin  found  that  the  range  of  current 
where  phase-slip  oscillations  occur  increases  as  ]  increases  and  for 
j  >  Jmax>  a  periodic  array  of  phase-slip  centers  can  exists.  Thus, 
the  dissipative  states  observed  experimentally  corresponds  to  the 
appearance  of  successive  PSCs.  They  found  that  the  linear  portion  of 
the  I-V  is  equivalent  to  that  of  a  segment  of  filament  of 
length  A  through  which  a  time-averaged  normal  current  of  magnitude 
j-  (3 Jmax  flows.  The  magnitude  of  the  excess  current  depends  on  ^  but 
is  independent  of  the  spatial  period  of  the  periodic  array  of  PSCs. 
e.g.  (1  =0.46  for  ^  =20,  (}=0.6  for  #  =10)-  On  the  other 

hand,  A  ~  7.5  £  (T)  independent  of  the  spatial  period  of  the  periodic 
array  and  independent  of  .  Figure  II. 3  shows  a  set  of  I-V  curves 
calculated  by  Kramer  and  Watts-Tobin. 

Deriving  the  theoretical  existence  of  an  array  of  PSCs  by  means  of 
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for  #  =10  and  #=20.  The  shaded  region  indicates  the  range  of  current  for  which 
Kramer  and  Baratoff  (reference  27)  found  phase-slip  oscillations,  d  is  the 
spatial  period  of  the  periodic  array  of  PSCs. 


a  theory  such  as  TDGL  is  a  big  advance  in  the  understanding  of  the 
features  in  the  I-V  curves  of  ODSFs.  However,  the  Kraner  and 
Watts-Tobin  equations  are  only  valid  if  £  (T)")^  (D  that  is,  in 

the  limit  of  local  equilibrium  or  slowly  varying  distribution.  In  tin, 
for  example,  this  limit  is  true  at  1-T/TC  <£<  0.003,  but  the  range  of 
temperature  in  which  voltage  steps  are  observed  is  much  wider. 
Furthermore,  this  model  was  developed  for  infinite  samples,  so  the  I-V 
curves  are  determined  by  the  separation  between  PSCs  and  not  by  the 
total  number  of  PSCs  activated. 
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One  important  qualitative  result  of  Kraner’s  work  ’  ’  ’  is  that 

during  the  phase-slip  oscillations,  the  magnitude  of  the  order  parameter 
is  depressed  during  the  whole  cycle  and  the  time  average  of  l IjJ i  is  much 
smaller  than  the  equilibrium  value.  This  result,  which  was  not 
anticipated  in  the  SBT  model,  has  been  observed  experimentally  as  will 
be  mentioned  in  section  II. 5. 

Octavio  and  Skocpol^0  solved  numerically  the  Kramer  and  Watts- 
Tobin  equation  for  a  finite  sample  .modified  by  inclusion  of  a  damping 
term  proportional  to  the  time  derivative  of  the  nonequilibrium 
potential.  They  also  found  that  the  order  parameter  is  depressed 
throughout  the  cycle,  and  moreover  that  the  oscillations  of  J  ,  Jn  and 
yKn-  at  the  center  of  the  PSC  propagate  like  damped  and  dispersive 
waves,  just  as  in  the  KSS  model1 

Kulik^1  has  developed  an  elegant  model  for  the  penetration  of  an 
electric  field  in  a  superconductor.  His  model,  like  the  SBT  model,  is 
based  in  the  dynamic  equations  of  the  two-fluid  model,  namely 
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C*'  C  n  =  0 

div  Pn  ♦  divps  =0 


(11.18) 


and  the  gauge  invariant  formulation  of  superconductivity: 

, 


(11.19) 


u  =  1  (ti  lip  +  2e4) 
1  2  dt 


where  yM  is  the  invariant  potential  of  the  superconductor  and  is 
equivalent  to  yM-3  in  the  SBT  model.  (|>  is  the  scalar  potential  of 

the  electric  field.  Kulik  calculates  the  electron  distribution  function 
using  the  Boltzmann  kinetic  equation,  neglecting  the  effect  of  the  gap, 
but  taking  into  account  the  scattering  of  the  normal  excitations  by  the 
collective  oscillations  of  the  order  parameter.  Kulik  gives  the  name 
"bogolons"  to  those  oscillations.  He  arrives  to  the  equation  for  the 
field  in  the  superconductor 


V 


(11.20) 


where  the  characteristic  length  J\  is  frequency  dependent  and  for  zero 
frequency  is  given  by 


A(oo=0)  =  (- 


1/2 


®  Nsviv, 


(11.21) 


where  n  is  the  total  density  of  electrons,  Ns  is  the  effective  density 
of  states  for  pair  electrons  and  is  approximately  equal  to  TI  A  J!(0)/2kT 
near  T  ,  N(0)  is  the  density  of  states  at  the  Fermi  surface,  1/VL  is 
the  elastic  (electron- impurity)  scattering  time,  1/V  is  the  relaxation 


23 


time  of  the  excess  charge  of  quasiparticles  and  is  determined  by  the 
electron-phonon  relaxation  time  '&£=  1/\)  ^  and  the  electron-bogolon 
relaxation  time  1/VS;  Vg=  l^g+Vp^. 

Kulik  found  that  JY(<jO=0)  increases  as  T  approaches  Tq  according 
to  -A(uj=0)  oC  (1-T/T  up  to  a  temperature  T.  at  which  V„ 

C  IS 

becomes  the  dominant  term.  _A.(u>=0)  is  constant  between  T,  and  T  . 

i  c 

According  to  Kulik,  not  only  tne  dynamic  oscillations  of  the  order 
parameter  but  also  any  static  inhomogeneities  of  the  order  parameter  in 
the  sample  will  contribute  to  the  conversion  of  the  charge  by  a  process 
similar  to  an  Andreev  reflection. 

For  a  one  dimensional  filament,  under  the  conditions  such  that  a 
PSC  is  present,  Kulik32  found  that  the  time-averaged  electric  field  E(x) 
that  determines  the  phase  jumps  decays  exponentially  from  the  center  of 
the  PSC  with  a  characteristic  length  _A(w  =  0)  given  by  equation 
(11.21).  The  total  potential  drop  across  the  phase-slip  center  is 
related  to  the  phase  jump  repetition  period  by  the  Josephson 
relation 


Kul ik32 


V  * 

defines  the  excess  current,  like  in  the  KSS  model 

* • 


(11.22) 
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as 


J  c  =  7M  J(e  (t))dt  (11.23) 

Ho 

where  0  is  the  phase  difference  across  the  filament.  However,  in  order 
to  determined  the  numerical  value  of  Jexc,  Kulik  takes  the  average  of 
the  relation  (II. 8),  like  in  the  SBT  model,  except  that,  contrary  to 


SBT,  he  takes  the  average  between  q=0  (Jg=0)  and  q=1/£  f3  ^Js=Jc^  and 
he  obtains  Jexc=5Jc/8.  This  result  appears  to  be  inappropriate  since 
the  average  of  J  should  be  taken  over  the  period  of  the  oscillation,  as 
indicated  by  equation  (11.23)  ,  i.e.  between  two  phase-slip  events  when 


For  a  periodic  array  of  n  PSCs  in  a  filament  of  length  L,  the  I-V 
curve  calculated  by  Kulik^  coincides  with  the  result  obtained 

qji 

previously  by  Tinkhair  for  the  generalization  of  the  SBT  model  for  n 


PSCs,  namely 


JW  =  Jexc  ♦ 


2nA{> 


'(—) 
\2n  JAS 


(11.24) 


Figure  II. 4  shows  the  theoretical  I-V  curves  obtained  by  Tinkhanr  . 

To  complete  this  review,  I  will  exanine  some  of  the  ideas  in  the 
Ivlev-Kopnin  model  of  a  PSC. 

According  to  Ivlev  and  Kopnin  ,  a  PSC  can  be  regarded  as  a 
topological  singularity  of  the  vortex  type  in  the  two  dimensional 
|x,ctj  space  where  they  form  a  periodic  structure.  The  electric  field 
satisfies  the  relation 


c  Edxdt  = 

Js  e 


(11.25) 


where  S  is  the  area  of  a  unit  cell  in  the  j  x,ctj  space.  Equation 
(11.25)  is  equivalent  to  the  Josephson  relation.  If  xq  is  the  length 
associated  to  the  PSC,  and  'J;  the  time  period  of  the  phase-  slip 
process,  then  the  average  voltage  across  one  PSC  is  V=  <E>  xQ  where 
<E)is  the  space  and  time  average  of  the  electric  field 


0  0.4  0.8  1.2  1.6  2.0  2.4 

x/le 

Figure  II. 4. Theoretical  I-V  curves  calculated  by  Tinkham  (reference  34) 
for  filaments  of  various  lengths  L  and  for  Iav„=0.5I„. 
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<E>  = 


1 


xo^ 


E  dx  dt  = 


V  = 


-'s 

nn-fi  nfito 


n  H  h 
exo\ 


V 


2e 


(11.26) 


(II. 27) 
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As  in  the  SBT  model,  Ivlev  et  al.  considered  the  structure  of  a 
PSC  in  two  length  scales,  one  region  smaller  than  (  where  the 

order  parameter  and  the  supercurrent  oscillate  in  time  and  another 
region  of  the  size  of  the  penetration  length  of  the  electric  field  X  ^ 
which  determines  the  relaxation  of  yks.  They  solved  numerically 

the  TDGL  equations  for  an  homogeneous  infinite  filament  like  in  the 
Kramer-  Baratoff  model,  but  in  the  limit  £  <(< 


Ivlev  et  al. 37,38 


found  periodic  oscillations  of  the  supercurrent 
1/2 


confined  to  a  region  0  <  x<  (  ,  but  they  found  that  the  region 

where  the  order  parameter  oscillates  is  much  smaller  than  f  (T). 

In  the  limit  5  £  ,  the  effect  of  the  alternating  supercurrent  on 

the  properties  of  the  PSC  can  be  disregarded  like  in  the  static 

pr»  pp  p£ 

approximation  of  Fink  y’  and  Galaiko  ■3”  °  and  the  chemical  potential 
of  the  pairs  to  each  side  of  the  PSC  can  be  taken  as  constant.  This 
static  approximation  is  valid  for  x  >>  (  '  and  there,  the  TDGL 

equations  become  a  system  of  time  independent  equations.  The  I-V  curves 


obtained  from  those  static  equations  start  with  an  infinite  slope  (for 

■31) 

an  infinite  sample)  at  j=j  ,  just  as  in  Tinkham's  model  . 

For  Jc  <"<  j  «  (  S  g-/  f  )^j  ,  the  I-V  curve  has  a  constant  slope  equal 
to  the  normal  resistance  and  stows  an  excess  current  independent  of  the 


voltage  at  T  close  to  T,,  jgxc=0.68  j  ,  but  it  depends  on  the  voltage  at 
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lower  temperatures. 

In  a  region  within  a  distance  from  the  center,  the  time 

average  of  the  modulus  of  the  order  parameter  is  much  smaller  than  the 
equilibrium  value,  just  as  in  Kramer's  model.  T  is  the  pair  breaking 
term  and  is  taken  to  be  F  C<  1  when  there  is  a  gap  in  the  energy 
spectrum,  which  is  the  case  in  the  Ivlev-Kopnin  calculations.  In  the 
region  x<<(£  ig.)^  the  oscillations  of  the  order  parameter 

and  of  the  nonequilibrium  potential  are  small,  but  the  oscillations  of 
the  supercurrent  are  still  considerable.  At  x  »  (  J  g.^-  )  ^  all 
quantities  become  time-independent.  Figure  II. 5  shows  the  approximate 
analytic  solution  of  the  dynamic  equations  corresponding  to  the  phase 
slippage  process  obtained  by  Ivlev  and  Kopnin^.  Those  solutions  are 
valid  in  the  temperature  range  where  the  phase-slip  frequency  CO  <^<  1/t>E 
and  the  pair  breaking  term  i . 


In  tin,  this  range  corresponds  to  0.0002  «  1-T/TC<<  0.005.  In  the 

limit  r«i  ,  the  penetration  length  of  the  electric  field  is,  according 

to  Ivlev  et  al.,  equal  to  O.^^f/P^. 

qq 

Ivlev  and  Kopnin  have  also  attacked  the  problem  of  calculating 

the  I-V  curve  of  a  ODSF  under  the  effects  of  microwave  irradiation.  In 

1  /2 

their  approximation,  they  assime  that  the  region  x<£  T  does  not 
contribute  to  the  calculation  of  the  I-V  curves  .  They  found  current 
steps  at  voltages  such  that 


£1=  muu 


(11.28) 


where  A  is  the  frequency  of  the  microwave  field,  ra  is  an  integer 
and  is  the  Josephson  frequency  of  the  phase-slip  process.  In  order 
to  calculate  the  magnitude  of  the  current  steps  it  is  necessary  to  know 
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the  exact  form  of  the  time  dependent  supercurrent.  Ivlev  and  KopninJ 
have  simplified  the  problem  by  assuming  that  the  space  and  time 
dependence  of  the  order  parameter  can  be  approximated  by 


A(x,t)0C 


tanh  (IXI/5P) 
tanh  (  x2(t)/$f2) 


if  |x|  >  x2(t) 
if  |x|  <  x2(t) 


(11.29) 


1  /O 

where  x2(t)  oscillates  from  zero  to  £  P  so  that 

tanh  x2(t)/£'f2  ~  x2(t)/£  J~2  .  They  found  that  the  size  of  the  steps  is 

given  by 


j 


=  2 


SrAfd 


A 


m 


(II. 30) 


CSnis  the  normal  conductivity  of  the  sanple,  E  is  the  d.c.  electric 
field  across  the  PSC,  is  the  differential  conductivity  dj/dE. 
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jr 
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e2P  E  E~  Ifml 

(Tc-T)1/2(u3m)5/2  r 


(II. 31) 


E«,  is  the  amplitude  of  the  a.c.  electric  field  in  the 
superconductor,  D  is  the  diffusion  constant  and  am  is  a  factor  that 
depends  on  the  form  of  X^t) 
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(x2(t’Vf/2)(1  -  cos((jumt'))dt’ 


(11.32) 
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These  results  are  a  very  rough  approximation  but  they  constitute 
the  first  attempt  to  explain  quantitatively  the  magnitude  of  the 
Josephson  steps  in  PSCs  taking  into  account  the  spatial  distribution  of 
the  a.c.  supercurrent 

As  we  have  seen  in  this  review,  the  problem  of  the  structure  of  a 

phase-slip  center  has  been  approached  in  different  ways.  Nevertheless, 

all  the  models  coincide  in  the  fact  that  the  structure  and  dynamics  of  a 

PSC  has  to  be  considered  on  at  least  two  length  scales:  a  short  region 

where  all  the  rapid  oscillatory  phenomena  occur,  and  a  longer  region 

determined  by  the  relaxation  of  the  nonequilibriun  potential.  In  the 

rest  of  this  report,  I  will  refer  to  those  regions  as  the  "core"  and  the 

1?  ?? 

"asymptotic  wings"  of  the  PSC.  Recently,  Baratoff  ’  introduced  a 
third  region  that  he  calls  the , "matching  region"  between  the  core  and 
the  wings  and  where  the  decay  of  the  nonequilibrium  potential  is 
non-exponential.  However,  since  no  observations  relating  directly  to 
this  region  have  been  made,  it  will  be  referred  to  only  infrequently  in 
the  rest  of  this  report.  In  fact,  since  the  simple  SBT  model  works  so 
well,  and  anticipates  most  of  the  general  features  of  the  subsequent 
more  microscopic  theories,  it  will  be  the  primary  descriptive  model  that 
is  used. 

Next, to  complete  this  review  of  earlier  work,  I  will  describe 
briefly  the  most  prominent  experiments  carried  out  on  PSCs  since  they 
were  first  observed  15  years  ago  until  the  present. 


II. 5  PREVIOUS  EXPERIMENTS  ON  PHASE-SLIP  CENTERS 


The  step-like  structure  connected  to  the  appearance  of  PSCs  was 
observed  for  the  first  time  in  the  R  versus  T  curves  of  tin  whiskers1. 

Later,  they  were  observed  in  the  I-V  curves  of  tin  whiskers  (1972), 

•a  m  c 

microbridges  of  tinJ  (1974),  aluminum  (1976),  indiun  (1976)  and 

niobium0  (1979)  and  also  in  the  I-V  curves  of  indium  whiskers  (1980). 

One  of  the  most  important  properties  of  those  I-V  curves  is  that  it 

is  possible,  in  principle,  to  determine  the  diffusion  length  of  the 

quasiparticle  charge  (or  characteristic  decay  length  of  the 

nonequilibrium  potential)  _A.q*  by  measuring  the  incremental 

differential  resistance.  According  to  the  SBT  model,  for  the  nth  step, 

we  have 


(11.33) 


where  wd  is  the  cross-sectional  area  of  the  filament  and  p  is  the 
normal  resistivity  of  the  sample. 

Many  experiments  on  PSCs  have  been  focused  primarly  in  determining 

_A.q„  by  means  of  the  equation (II. 33).  However,  the  values  obtained  by 

using  this  method  differ  from  the  true  values  because  the  heat  generated 

in  the  PSCs  adds  some  curvature  to  the  I-V  curves  and  the  measured 

values  of  dV/dl  are  overestimated.  This  was  the  case  of  the 

experimental  results  on  tin  microbridges  obtained  by  SBT^  who  were 

unable  to  detect  the  divergence  at  Tc  of  Aq*  and  consequently,  they 

concluded  that  the  inelastic  scattering  time  lj>  c  rather  than  the  branch 

mixing  time  determined  the  appropriate  diffusion  length  of  the 
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quasiparticles.  It  had  been  shown  theoretically  ’ 


that  the 
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relaxation  time  of  the  quasiparticle  charge  imbalance  mode  of 
disequilibrium  'Sq*,  which  determines  -Aq*,  diverges  at  Tc  as 
oC(1-T/Tc)"1/2,  so  that  A q t  should  diverge  as  A^oCtl-T/^)'1^. 

This  contradiction  between  experimental  results  and  theoretical 
prediction  puzzled  many  researchers  for  some  time.  The  first 
observation  of  the  divergence  at  Tc‘of  -Aq*  was  reported  by  Jillie^  who 
measured  the  differential  resistance  of  the  I-V  curves  of  indium 
microbridges.  This  result  raised  the  question  of  whether  this  effect 

was  dependent  on  the  superconductor  used. 

40 

Kadin  et  al.  subsequently  showed  how  heating  effects  can  obscure 
the  divergence  of  the  differential  resistance  as  T  approaches  T  ,  and 
after  correcting  their  data  for  heating  effects  they  were  able  to 
observe  the  temperature  and  magnetic  field  dependence  of  the 
differential  resistance  of  tin  microbridges  in  the  presence  of  a 
magnetic  field.  Their  results  agree  well  with  the  transverse  mode 
relaxation  time  in  the  presence  of  a  pair-breaking  mechanism  as  derived 
by  Schmid  and  Schon1^. 

It  is  also  possible  to  measure  directly  the  potential  of  the  pairs 
and  the  potential  of  the  quasiparticles  by  attaching  sets  of 
tunnel  junction  probes  on  a  microbridge.  The  first  of  this  type  of 
measurements  were  carried  out  by  Dolan  and  Jaokel11.  They  attached  a 
set  of  silver  probes  and  a  set  of  lead  probes  to  a  long  tin  microbridge. 
They  measured  the  voltage  at  different  probes  and  found  that  the 
nonequilibrium  potential  decayed  exponentially  away  from  the 

center  of  the  PSC,  supporting  the  description  of  the  PSC  according  to 
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the  SST  model.  Dolan  and  Jackel  also  observed  the  expected  temperature 

divergence  of  this  directly  measured  decay  length.  Their  data 

corresponds  to  a  relaxation  time  ^Q#=(0.9i0.2)x10"10sec(1-T/Tc)’V2. 

Using  a  similar  geometry,  Skocpol  and  Jackel  observed  the 

reduction  of  the  energy  gap  at  the  core  of  the  PSC.  Their  samples  were 

tin  microbridges  with  lead  tunnel  junction  probes.  By  measuring  the  I-V 

curves  of  the  tin/tin  oxide/lead  junctions,  they  were  able  to  obtain  the 

value  of  the  energy  gap  along  the  microbridge.  They  found  that  the 

first  PSC  appears  where  the  gap  is  minimum  and,  once  the  PSC  appears, 

the  gap  is  depressed  by  more  than  90  %.  This  result  confirmed  the 

models  of  Kramer  et  al.2^”2^  Ivlev  et  al.^“^  and  Octavio  et  al.^°. 
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Figure  II. 6  shows  the  experimental  results  of  Skocpol  and  Jackel  . 

42 

Stuivinga  et  al.  prepared  long  aluminum  strips  with  normal  copper 

probes,  and  they  measured  the  decay  with  distance  from  the  center  of  the 

PSC  of  the  normal  potential.  They  found  that  the  relaxation  time 

-1/2 

diverges  at  Tc  as  (1-T/TC)  .  This  is  the  expected  temperature 

dependence  if  inelastic  (electron-phonon)  scattering  time  is  the 
predominant  mechanism  in  relaxing  the  charge  imbalance.  However, 

llg 

Stuivinga  et  al.  had  predicted  that  in  materials  like  Al  and  Zn,  where 

the  inelastic  scattering  time  'Sg.  is  relatively  long,  the  relaxation  of 

the  charge  imbalance  is  dominated  by  the  pair  breaking  effect  of  the 

3upercurrent  and  therefore,  the  measured  relaxation  time  should  diverge 

faster  than  (1-T/T  )  as  T  approaches  T  .  They  explained  this 
c  c 

discrepancy  between  their  theoretical  model  and  their  experimental 
results  in  terms  of  the  dependence  with  position  (because  of  Js(x))  of 


3' 


Figure  I I. 6. Experimental  measure 
microbridge  obtained 
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the  relaxation  time  associated  with  that  mechanism. 

The  oscillations  of  the  order  parameter,  or  more  precisely,  the 
oscillations  of  the  supercurrent  in  a  PSC  have  been  detected  in 

■3 

different  ways.  Skocpol  et  al.  (SBT)  observed  the  current  steps  in 

the  I-V  curves  of  long  (120 m)  tin  microbridges  irradiated  with 

microwaves  (f  =  10  GHz)  at  voltages  corresponding  to  n  times 

(n=1 ,1/2, 1/3)  the  Josephson  voltage  in  the  first  and  second  PSC. 

411 

Tidecks  and  Minnigerode  observed  Josephson  steps  in  the  I-V  curves  of 
tin  whiskers  (500  long)  irradiated  with  microwaves  of  various 
frequencies  (from  f=0.2  GHz  to  f=1  GHz).  They  proposed  an  expression 
for  the  magnitude  of  the  steps  using  an  approximation  for  the  spatially 
averaged  time  dependent  supercurrent 

rs(t)=  ^IC(1  +  c°s<p(t))  (11.34) 

(p( t)=  <p(0)  +  V(t’)dt'  (11.35) 


Their  result  for  the  time  average  of  the  supercurrent  is 

<Is(t)>  =-Lyi_2e_  y aC0S(cp  (O)-(p(a) ) 

2  2hu> 


(11.36) 


where  Vg  and  ( pg  are  the  amplitude  and  phase  of  the  microwave  voltage  in 
the  sample. 

According  to  equation  (11.36),  the  magnitude  of  the  steps, 
normalized  to  the  critical  current,  is  proportional  to  the  microwave 
voltage.  However,  no  comparison  between  this  model  and  the  experimental 


results  is  discussed  in  reference  44. 
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Josephson  steps  have  also  been  observed  in  the  I-V  curves  of  "long" 

Nb  nanobridges  (from  L=0.12  f<  m  to  L=0.9  fk  m)  by  Laibowitz  et  al.^  and 

7 

also  in  the  I-V  curves  of  In  whiskers  by  Tiedecks  et  al.  . 

41 

In  a  different  type  of  measurement  ,  Skocpol  and  Jackel  have 
detected  the  oscillations  of  the  pair  potential  jk^  in  a  PS C  by 
observing  small  current  steps  induced  in  the  I-V  curves  of 

superconducting  tunnel  probes  located  in  the  vicinity  of  a  PSC  in  a  tin 

microbridge. 

Other  experimental  evidences  that  prove  the  oscillatory  nature  of 

PSCs  have  been  the  observation  of  voltage  locking  between  two  PSCs  and 

the  detection  of  standing  charge  imbalance  waves.  The  latter  effect  was 

uc 

observed  by  Kadin  et  al.  in  In  microbridges  of  various  lengths  and  it 

occurs  when  the  wavelength  of  the  charge  imbalance  wave  "X  and  the 

length  of  the  microbridge  L  are  such  that  2L=(n+1/2)\  where  n=0,1,... 

Voltage  locking  occurs  when  the  waves  emitted  by  one  PSC  modulate 
the  voltage  driving  the  other  PSC.  This  coherent  behavior  has  been 
detected  by  Lindelof  et  al.^’^  who  measured  the  radiation  emitted  by 
pairs  of  short  microbridges  made  of  tin  or  indium  with  various 

48  49 

separations  between  the  microbridges,  and  also  by  Jillie  et  al.  ’  who 
measured  directly  the  I-V  curves  of  pairs  of  short  indiun  microbridges. 

According  to  SET,  PSCs  occur  where  the  superc^-rent  I  reaches  the 
critical  value  I0(x)  at  some  point  in  the  filament.  Vfoen  a  PSC  is 
already  present,  the  supercurrent  is  depressed  in  its  neighborhood, 
I3(x)=  I-I  (x)  (see  figure  II. 1).  In  consenquence,  successive  PSCs  are 
more  likely  to  occur  far  from  the  ones  already  present.  This  apparent 
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enhancement  of  the  critical  current  in  the  vicinity  of  a  PSC  was 

50 

observed  by  Meyer  et  al.  in  the  I -V  curves  of  tin  whiskers  with 
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potential  probes  and  also  by  Lindelof  and  Hansen  in  pairs  of  tin  and 
indiun  short  microbridges. 

This  interaction  via  d.c.  quasiparticle  current  was  studied  in 

Q 

more  detail  by  Jillie  et  al.  in  pairs  of  short  indium  microbridges. 

They  measured  the  change  of  the  apparent  critical  current  of  one  bridge 
(detector)  as  a  function  of  the  current  and  voltage  in  the  other  bridge 
(source)  in  the  cases  of  both  opposing  and  series  current  bias.  When 
the  current  flows  in  the  same  direction  in  both  microbridges,  an 
increase  in  the  apparent  critical  current  of  the  detector  should  be 
observed  according  to  SET.  When  the  current  flows  in  opposite 
directions,  the  supercurrent  in  the  detector  is  increased  and  the 
apparent  critical  current  of  the  detector  decreases.  This  change  in  the 
apparent  critical  current  is  equal,  in  magnitude,  to  the  quasiparticle 
current  generated  in  the  source  that  flows  in  the  detector.  On  the 
other  hand,  the  Joule  heat  generated  in  the  source  decreases  the 
critical  current  of  the  detector  independently  of  the  relative  direction 
of  the  currents. 

q 

Jillie  et  al.  separated  those  two  effects  by  decomposing  the 
change  of  the  critical  current  of  the  detector  in  parts  which  are 
symmetric  and  antisymmetric  respect  to  the  relative  direction  of  the 
currents.  They  found  that  these  effects  vanished  for  separations 
greater  than  4  fKm.  This  happens  because  in  their  geometry,  the 
two-dimensional  spreading  of  the  nonequilibrium  quasiparticle  current 
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weakens  the  effect.  They  did  not  observe  an  increase  in  the  apparent 
critical  current  of  the  detector  because  of  an  overwhelming  heating 
effect.  From  their  data,  they  estimated  a  temperature  rise  in  the 
detector  of  1.3  mK  per  each  nW  of  power  dissipated  in  the  source.  In 
the  work  to  be  reported  here,  these  difficulties  were  overcome  by  use  of 
a  more  favorable  sample  geometry. 


CHAPTER  III 


EXPERIMENTAL  TECHNIQUES 

III. 1  FABRICATION  OF  SAMPLES 

Three  different  types  of  samples  were  prepared  during  the  course  of 
this  work.  I  will  refer  to  them  as  samples  of  type  M,  samples  of  type 
S,  and  samples  of  type  T. 

Samples  of  type  M,  used  to  study  the  effects  of  an  external 
microwave  radiation  on  a  PSC,  are  tin  microbridges  of  various  lengths, 
from  l^is  to  UOy^m,  and  about  l^m  wide.  Figure  III.1  shows  four 
samples  of  type  M.  Samples  of  type  S,  used  to  study  the  interaction 
between  two  PSCs,  consist  of  a  tin  microbridge  40^  m  long  and  about  lyum 
wide  with  four  side  probes,  all  part  of  the  same  film.  The  distance 
between  the  probes  was  chosen  to  be  of  the  order  of  the  quasiparticle 
diffusion  length  so  that  only  one  PSC  nucleates  in  any  of  the  five 
segments  of  the  bridge  when  current  is  fed  through  the  adjoining  probes. 
The  probes  were  made  wider  than  the  bridge  to  allow  application  of 
currents  well  above  the  critical  current  of  the  bridge  without  loss  of 
superconductivity  in  the  probes.  Figure  III. 2  shows  two  samples  of  type 
S.  Samples  of  type  T,  used  to  measured  the  nonequilibriun  potential  in 
a  PSC,  also  consist  of  a  tin  microbridge  40^.m  long  and  about  I^Mm  wide, 
but  rather  than  side  probes,  they  have  four  tunnel  junction  probes.  In 
some  of  the  samples  of  type  T,  the  middle  segment  of  the  microbridge  was 
made  narrower  than  the  rest  of  the  microbridge  in  order  to  facilitate 
the  nucleation  of  a  PSC  in  that  location.  This  turned  out  to  be  a 
successful  technique,  as  will  be  discussed  later.  Figure  III. 3  shows 
two  samples  of  type  T. 
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Figure  III. 1. SEW  pictures  of  samples  (a)  M-29,  (b)  M-15,  (c)  M-30  and 
(d)  M-25.  Sample  M-25  was  made  with  the  same  mask  used  to 
make  samples  of  type  T,  but  doing  a  single  tin  evaporation. 
The  microbridge  is  the  horizontal  strip. 


V, 
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All  our  samples  were  fabricated  using  photolithographic  techniques 
which  I  will  describe  in  the  following  sections,  beginning  with  the 
substrate  cleaning  procedures. 

III. 1.1  SUBSTRATE  CLEANING  PROCEDURE 

Anyone  involved  in  microfabrication  and  thin  film  technology  knows 
that  there  is  no  universal  recipe  for  cleaning  substrates.  It  depends 
on  the  substrate  material,  the  purpose  of  the  cleaning  and  the  time  that 
one  is  willing  to  spend  cleaning.  We  have  combined  several  existing 
methods  to  make  a  simpler  one  that  we  found  to  be  effective  for  cleaning 
our  glass  substrates. 

Our  substrates  were  made  from  microscope  slides  1  mm  thick  'which 
had  been  cut  into  pieces  of  dimensions  1"x1/4",  1"x1/2"  or  1"x1";  this 
means  that  our  substrates  are  subjected  to  considerable  handling  before 
they  are  cleaned.  We  start  oy  wiping  the  substrate  with  a  piece  of 
lint-free  wiping  cloth  wet  with  ethanol  in  order  to  remove  the  grossest 
dirt.  We  used  ethanol  because  it  did  not  dissolve  the  disposable 
gloves . 

The  next  step  is  to  boil  the  substrate  in  trichloroethylene  (TCE) 
for  five  minutes.  TCE  dissolves  most  of  the  organic  contaminants.  This 
step  is  followed  by  sequential  ultrasonic  agitation  in  TCE,  acetone  and 
methanol  for  approximately  five  minutes  in  each  solvent.  Finally,  the 
substrates  are  dried  with  filtered  pure  nitrogen  and  then  they  are 
stored  in  tightly  sealed  containers. 

As  will  be  discussed  in  section  III. 1.3,  during  the  fabrication  of 
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our  masks  we  need  to  prepare  chrome-coated  1"xV’  glass  slides.  We  found 
that  in  order  to  improve  the  quality  of  the  Cr  films,  in  particular  to 
avoid  pinholes,  it  was  necessary  to  include  an  additional  step  in  the 
cleaning  process.  Before  the  substrates  were  boiled  in  TCE,  they  had  to 
be  boiled  in  a  solution  of  about  1  part  of  Micro  in  10  parts  of  water 
for  about  30  minutes  and  then  rinsed  thoroughly  in  distilled  water  and 
dried  with  pure  nitrogen;  the  other  steps  were  then  followed.  (Micro 
is  a  liquid  laboratory  cleaner  manufactured  by  International  Products 
Corp.) . 

III. 1.2  PREPARATION  OF  THE  PHOTORESIST  STENCIL 

In  order  to  prepare  our  samples  we  used  a  lift-off  technique 
SI  B2 

developed  by  Dunkleberger  and  Dolan  .  Their  method  consists  of 
fabricating  two  layers  of  photoresist  separated  by  a  thin  film  of 
aluminum.  The  pattern  is  defined  in  the  top  layer  of  photoresist,  which 
is  offset  from  the  substrate  by  the  bottom  layer  of  photoresist.  We 
used  Shipley  AZ1350J  photoresist  to  prepare  the  stencils. 

The  clean  substrate  is  placed  in  a  photoresist  spinner,  covered 
with  a  few  drops  of  photoresist  and  spun  at  5500  RPM  for  25  seconds  to 
obtain  a  1,5^m  thick  layer  of  photoresist.  The  substrates  covered  with 
photoresist  are  then  baked  at  95°Ci5*C  for  25  minutes  and  exposed  with 
ultraviolet  light  for  about  3  minutes  to  make  the  photoresist  completely 
soluble  in  the  developer.  We  usually  wait  more  than  10  minutes  between 
the  bake  and  the  exposure  in  order  to  allow  the  photoresist  to  recover 
the  moisture  lost  during  the  bake. 
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A  800  X  to  1000  A  film  of  aluminun  is  deposited  on  top  of  the 
exposed  photoresist  and  then  a  second  layer  of  1.5y4m  of  photoresist  is 
spun  on  top  of  the  A1  film.  This  second  layer  has  to  be  baked  at  a 
tsnperature  a  few  degrees  lower  than  the  temperature  at  which  the  first 
layer  is  baked  (which  is  higher  than  the  90  °C  temperature  recommended  by 
Shipley),  otherwise,  lumps  on  the  surface  of  the  aluminun  film  appear 
due  to  the  release  of  gases  from  the  bottom  layer  of  photoresist  during 
the  second  bake.  The  sandwich  (exposed  p'notoresist)/aluminum/(unexposed 
photoresist)  is  now  ready  to  be  used. 

Contact  printing  was  used  to  copy  the  pattern  from  the  mask  onto 
the  top  layer  of  photoresist.  The  masks  consist  of  a  chrome  film 
patterned  on  a  2"x2"  piece  of  conformable  glass  0.25  mm  thick.  The 
substrate  coated  with  photoresist  is  placed  in  a  hole  in  a  plexiglass 
block.  The  chrome  mask  is  held  by  suction  in  intimate  contact  with  the 
substrate.  The  exposure  time  was  90  seconds  for  sanples  placed  2  feet 
away  from  the  ultraviolet  (UV)  light  lamp.  The  portions  of  photoresist 
exposed  to  the  UV  light  are  then  removed  by  soaking  the  sample  in  the 
developer  (1  part  Shipley  Microposit  351  in  5  parts  of  distilled  water). 
The  photoresist  is  developed  for  approximately  60  seconds. 

The  aluminum  film  is  etched  in  a  solution  of  84%  H^PO^,  5%  H20,  5% 
CH^COOH  and  5%  HNO^  by  volume.  The  bottom  layer  of  photoresist  is 
partially  removed  by  immersing  the  sample  in  the  developer  for  only  a 
few  seconds  (usually  between  10  and  13  seconds)  in  order  to  create  an 
undercut  under  the  top  layer.  Finally,  the  sample  is  immersed  in  the 
aluminum  etch  in  order  to  get  rid  of  the  aluminun  film  under  the 
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overhanging  edge  of  photoresist,  by  doing  so  the  undercut  can  be 
inspected  by  looking  from  the  top  using  an  optical  microscope.  Figure 
III. 4  shows  photoresist  stencils  ready  to  be  used  to  prepare  (a)  a 
sample  of  type  M  and  (b)  a  sample  of  type  T. 

Figure  111.5(b)  shows  a  close-up  of  the  stencil  shown  in  figure 
111.4(b).  The  undulations  in  the  surface  of  the  photoresist  in  figure 
111.5(a)  are  produced  by  standing  waves  formed  when  there  is  a 
separation  between  the  mask  and  the  surface  of  the  photoresist  during 
the  contact  printing.  Very  often,  residues  of  photoresist  on  the  chrome 
mask  prevent  a  good  contact,  for  that  reason  it  is  necessary  to  inspect 
the  mask  carefully  and  to  clean  it  in  case  that  it  is  not  perfectly 
clean  before  doing  the  contact  printing.  Good  contact  between  the  mask 
and  the  photoresist  was  found  to  be  more  critical  than  the  exposure  time 
and  the  developing  time  for  obtaining  vertical  walls  in  the  pattern  of 
photoresist. 

The  mask  used  to  make  the  samples  of  type  M  and  the  samples  of  type 
S  was  made  by  a  commercial  firm  (Advance  Reproductions  Corporation. 

North  Andover,  Massachusetts).  A  photograph  of  the  entire  mask  is  shewn 
in  figure  111.6(a).  The  dark  areas  are  the  areas  covered  with  chrome, 
the  white  areas  correspond  to  the  areas  that  become  metallized  during 
the  deposition  of  the  samples.  The  mask  consists  of  six  sections.  The 
widest  of  them  is  the  mask  used  to  make  the  samples  of  type  S  and  a 
blow-up  of  its  middle  part  is  shown  in  figure  111.6(d).  Each  one  of  the 
other  five  sections  contain  four  microbridges  of  various  lengths. 

Figures  111.6(b)  and  111.6(c)  show  a  blow-up  of  two  of  those  masks. 
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Figure  III.4.SEM  pictures  of  photoresist  stencils  ready  to  be  used  to 
prepare  (a)  a  sample  of  type  M  and  (b)  a  sample  of  type  T. 


Figure  III. 5. Close-up  of  two  stencils  used  to  prepare  samples  of  type  T. 

We  can  see  the  difference  between  a  bad  (a)  and  a  good  (b) 
contact  printing.  See  text  for  details.  Both  pictures  were 
taken  at  the  sane  magnification. 


50 


From  the  original  or  master  mask,  we  made  separate  copies  of  each  one  of 
those  six  sections  in  individual  2"x2"  plates  of  conformable  glass 
precoated  with  chrome. 

The  narrowest  linewidth  in  the  original  masks  is  2/4.  m  which  becomes 
wider  during  the  processes  of  making  the  copy  in  the  conformable  glass 
and  preparing  the  photoresist  stencil.  We  can  see  in  figure  III. 4(a) 
that  tre  line  has  become  3/4 m  wide.  This  widening  of  the  line  is 
compensated  by  doing  an  oblique  evaporation  for  the  deposition  of  the 
samples.  During  the  evaporation  the  substrate  is  held  at  an  angle 
relative  to  the  horizontal  plane.  By  doing  this,  the  width  w  of  the 
bridge  fabricated  is  considerably  narrower  than  the  mask  opening.  From 
figure  III. 7  we  can  see  that  w=a-htan 0  .  This  technique  also  allows  us 
to  produce  samples  of  type  S  in  which  the  bridge  is  made  narrower  than 
the  probes.  More  details  of  the  evaporation  techniques  are  given  in 
section  III. 1.4;  but  first,  I  would  like  to  explain  briefly  how  we  made 
the  mask  used  to  produce  the  samples  of  type  T. 

III. 1.3  FABRICATION  OF  MASKS 

As  it  was  mentioned  in  the  previous  section,  the  mask  used  to  make 
samples  of  type  S  and  samples  of  type  M  was  made  by  a  commercial  firm. 

At  that  time  we  did  not  have  the  expertise  nor  the  facilities  to  produce 
our  own  masks.  We  had  tried  making  masks  of  aluminim  film  and  using  a 
diamond  knife,  but  we  found  that  it  was  a  time-con sun  ing  process  with 
little  success  and  the  masks  produced  did  not  last  long  due  to  the 
softness  of  the  aluminum  film  and  their  poor  adhesion  to  the  glass 
substrates.  As  the  need  for  producing  samples  of  diverse  geometries  for 
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different  members  of  our  group  became  more  apparent,  and  we  found  that 

ns 

we  could  not  afford  to  order  new  commercial  masks,  W.  C.  Dan chi  set 
up  a  system  to  produce  our  own  masks.  With  a  light  board,  a  few  sheets 
of  rubylith,  a  few  sheets  of  clear  acetate,  black  tape  of  various  widths 
and  some  glass  plates  precoated  with  chrome  and  photoresist,  he  was  able 
to  produce  masks  with  lines  as  narrow  as  5^m  by  making  sucessive 
reduction  of  a  pattern. 

The  first  step  is  to  make  a  pattern  in  rubylith  or  with  black  tape. 

(Rubylith  is  a  material  consisting  of  a  thin  red  film  laminated  with  a 

heavier  clear  backing  so  that  the  red  film  can  be  scribed  and  peeled 

away,  leaving  the  desired  pattern).  This  pattern  is  reduced  25  times  by 

taking  a  picture  with  a  35  mm  camera.  The  35  mm  negative  is  used  to 

make  a  copy  by  contact  printing  on  a  1"x1"  glass  slide  coated  with 

chrome.  (Those  are  the  1"x1"  glass  slides  that  require  the  additional 

cleaning  step  mentioned  in  section  III.1).  The  next  step  consists  of 

doing  a  lOx  reduction  of  the  chrome  copy  using  an  optical  microscope. 
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This  second  step  was  possible  thanks  to  D.  J.  Frank  who  designed  a 
projection  attachment  to  use  with  our  Zeiss  microscope.  By  using 
different  objectives,  it  is  possible  to  do  lOx,  20x,  50x  and  lOOx 
reductions.  ^ 

i 

In  order  to  make  the  samples  of  type  T,  it  is  necessary  to  have  a 
ma3k  with  lines  2^111  wide  or  narrower,  which  means  that  we  had  to  use 
the  objective  with  the  highest  resolution.  On  the  other  hand,  the  area 
illuminated  gets  snaller  as  the  power  of  the  objective  is  increased. 

For  that  reason  we  had  to  make  our  masks  in  three  separate  reductions. 
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The  part  of  the  mask  that  defines  the  pads  is  made  by  reducing  the 
rubylith  pattern  25  times  using  the  35  mm  camera.  Then,  a  25x10x 
reduction  is  made  to  produce  the  intermediate  region  that  connects  the 
pads  to  the  smallest  region,  and  finally  a  25x50x  reduction  is  made  to 
produce  the  narrowest  lines.  We  obtained  lines  as  narrow  as  1.5yHm 
using  this  technique.  Figure  111.8(a)  shows  a  picture  of  one  of  those 
masks.  Figures  111.8(b)  and  111.8(c)  show  a  blow-up  of  the  central  part 
of  the  mask  for  two  different  masks,  both  used  to  make  samples  of  type 
T,  but  one  (b)  with  a  microbridge  of  uniform  width  and  the  other  (c) 
with  a  microbridge  with  a  narrow  segment  in  the  middle. 

III. 1.4  EVAPORATION  OF  THE  SAMPLES 

Most  of  our  tin  sanples  were  evaporated  onto  glass  substrates 
maintained  at  liquid  nitrogen  temperature,  and  a  few  samples  were 
deposited  at  roan  temperature.  We  did  not  find  any  difference  in  the 
quality  of  those  films  as  long  as  they  were  deposited  at  a  very  high 
rate,  200  A/sec  or  higher. 

The  vacuum  system  consists  of  a  mechanical  pump  with  a  molecular 
sieve  trap  and  a  cryopump,  with  an  ultimate  pressure  of  about  10~^  Torr. 
Some  of  our  earliest  samples  of  type  M  were  made  in  the  "old"  evaporator 
with  a  diffusion  pump  instead  of  a  cryopunp,  the  ultimate  pressure  being 

_7 

also  about  10  Torr. 

Samples  of  type  M  and  samples  of  type  S  were  made  in  a  single 
evaporation  of  a  tin  film  1000  A  thick  with  the  substrate  forming  an 
angle  9  with  the  horizontal  plane  as  indicated  in  figure  III. 7. 

O 

In  order  to  prepare  samples  of  type  T,  we  evaporated  a  1000  A  thick 
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Figure  III. 8. (a)  Micrograph  of  one  mask  used  to  make  samples  of  type  T. 
(b)(c)  Close-ups  of  the  central  part  of  the  masks. 
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tin  film,  oxidized  it  in  a  glow  discharge  of  pure  oxygen,  and  after 
rotating  the  substrate  to  the  correct  angle,  a  layer  of  normal  metal  was 
deposited.  The  glow  discharge  was  maintained  for  a  period  of  between  30 
seconds  and  100  seconds,  depending  on  the  desired  barrier  resistance. 

_p 

The  initial  pressure  of  the  oxygen  was  3x10  Torr.  We  used  silver  or 
indiun  as  the  normal  electrodes.  (The  indium  is  normal  in  the  range  of 
temperature  in  which  the  measurements  were  made) .  The  samples  are 
removed  from  the  evaporator  and  acetone  is  sprayed  on  them  to  dissolve 
the  unexposed  photoresist  and  lift  off  the  unwanted  metal. 

Next,  I  would  like  to  describe  other  methods  of  fabrication  that  we 
tried  during  the  progress  of  this  work. 

III. 1.5  OTHER  METHODS  OF  FABRICATION 

The  diamond  knife  technique  was  the  first  technique  that  we  used  to 
fabricate  microbridges.  In  order  to  make  cuts  on  a  thin  film,  the 
substrate  is  slid  beneath  a  diamond  knife  that  can  be  moved  vertically, 
so  that  it  is  possible,  in  principle,  to  control  the  pressure  with  which 
the  knife  is  placed  a£  \inst  the  film.  The  substrate  can  be  moved 
horizontally  in  two  perpendicular  directions  by  means  of  differential 
screws  attached  to  the  translation  stage,  so  that  parallel  cuts  can  be 
made,  like  those  shown  in  figure  111.9(b),  (c)  and  (d) ,  creating  a  long 
constriction  connecting  two  large  areas  of  the  film.  To  produce  short 
microbridges,  like  the  one  shown  in  figure  111.9(a),  movement  of  the 
substrate  along  only  one  direction  is  needed. 

As  we  can  see  in  figure  III. 9,  the  width  and  the  depth  of  the  cut 
is  different  for  different  samples,  this  happens  because  it  is  very 
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difficult  to  see  when  the  knife  is  touching  the  film  and  therefore,  it 
is  hard  to  decide  how  much  the  knife  can  be  moved  downward.  Sometimes 
the  position  of  the  knife  relative  to  the  film  changed  during  the 
dragging,  and  as  a  consequence,  the  film  was  not  cut  completely  leaving 
metallic  shorts  in  parallel  with  the  microbridge.  Other  inconveniences 
with  the  diamond  knife  technique  are:  firstly,  the  fact  that  the 
microbridges  produced  by  this  method  are  thicker  at  the  edges  because  of 
the  metal  pushed  by  the  knife  during  the  cuts,  and  secondly,  the  variety 
of  geometries  that  are  possible  to  fabricate  is  very  limited.  More 
details  about  the  diamond  knife  technique  of  microfabrication  can  be 
found  in  references  55  and  56. 

Another  technique  that  was  used  but  abandoned  later  consists  of  the 
lift-off  of  a  single  layer  of  photoresist  treated  with  chlorobenzene. 
This  technique  was  developed  by  Hatzakis  et  al.  '  and  it  consists  of 
soaking  the  baked  photoresist  in  chlorobenzene  either  before  or  after 
the  exposure.  (I  found  that  this  method  works  best  when  the  soak  is 
done  after  the  exposure) .  This  process  strengthens  the  surface  of  the 
photoresist  making  it  harder  to  dissolve  in  the  developer  and 
consequently,  an  undercut  profile  is  obtained.  Figure  III. 10  shows  SEM 
pictures  of  two  profiles  of  AZ1350J  that  have  been  soaked  in 
chlorobenzene  for  (a)  5  minutes  and  (b)  10  minutes  after  the  exposure. 
The  thickness  of  the  hardened  layer  of  photoresist  agrees  well  with  the 
results  reported  in  reference  57.  This  technique  works  very  well  to 
produce  clean  lift-offs;  however,  the  photoresist  stencils  can  not  be 
used  for  oblique  evaporations 
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Figure  III.10.SEM  picture  of  stencils  of  AZ135QJ  photoresist  treated 
with  chlorobenzene. 
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III. 2  MEASUREMENT  TECHNIQUES 

In  order  to  attach  leads  to  the  samples,  indiun  patches  were 
applied  to  the  substrates  before  the  evaporation.  After  the  samples 
were  made,  wires  were  connected  to  the  indiun  patches  using  a 
battery-powered  soldering  iron.  In  samples  of  type  S  and  sample  of  type 
T,  one  current  lead  and  one  voltage  lead  are  attached  to  each  end  of  the 
microbridge  and  to  each  one  of  the  four  probes. 

Measurements  on  samples  of  type  M  were  taken  using  a  two-terminal 
or  a  four-terminal  configuration,  depending  on  what  kind  of  set-up  was 
used  to  apply  microwaves  to  the  samples.  In  some  experiments,  the 
samples  were  placed  in  an  X-band  waveguide  at  a  distance  equal  to  A  /4 

o 

fhom  the  closed  end.  In  those  experiments  the  samples  were  two-probed 
because  the  slit  through  which  the  sample  is  inserted  in  the  waveguide 
is  almost  as  narrow  as  the  thickness  of  the  glass  substrates.  In  other 
experiments,  the  microwave  radiation  was  applied  to  the  sample  using  a 
semi-rigid  coaxial  cable  of  0.141"  outer  diameter  terminated  in  an 
inductive  loop.  This  loop  was  placed  parallel  to  the  sample  at  a 
distance  about  3  mm  from  it.  In  those  experiments  the  samples  were 
four-probed.  We  used  BNC  connectors  to  adapt  the  waveguide  to  the 
coaxial  cable  and  we  found  experimentally  that  the  microwave  radiation 
is  attenuated  by  20  dB  when  the  BNC-coaxial-loop  configuration  is  used 
instead  of  the  waveguide. 

For  all  the  samples,  the  voltage  leads  are  continuous  copper  wires 
connecting  the  samples  directly  to  a  Keithley  Nanovoltmeter.  The 
current  leads  are  connected  to  a  ten-pin  connector  and  from  there  to 
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battery-powered  current  sources.  Up  to  this  point  the  banks  have  been 
shorted  together  by  wide  segments  of  the  film  which  are  now  scraped  with 
a  Circon  microknife.  The  sanples  are  tested  for  continuity  at  room 
temperature  by  applying  a  small  current,  usually  less  than  5y<A,  and 
measuring  the  voltage  with  the  Nanovoltmeter. 

The  samples  are  immersed  directly  into  the  liquid  helium  bath  whose 
temperature  was  controlled  using  a  pressure  regulator.  During  the 
measurements,  the  temperature  was  kept  constant  within  one  millidegree 
for  up  to  five  hours,  as  confirmed  by  a  germanium  thermometer.  Figure 
III. 11  shows  a  schematic  diagram  of  the  circuit  used  to  measure  the  I-V 
curves  of  the  samples. 

Once  the  fabrication  of  the  samples  was  completed,  the  most  serious 
problem  encountered  was  that  of  burnout,  meaning  the  destruction  of  the 
microbridge  by  accidental  electric  pulses  caused  by  turning  electronic 
equipment  on  and  off,  or  by  simply  touching  the  leads  connected  to  the 
sample.  Sometimes,  samples  burned  out  without  apparent  reason.  Figure 
III.  12  shows  different  levels  of  destruction,  from  areas  about  1  {jk m)*-, 
to  enormous  catastrophes  where  the  area  destroyed  is  as  big  as 

4  2 

10  (^m)  .  The  results  presented  in  this  report  were  obtained  with  the 

samples  that  survived  for  the  longest  period  of  time,  and  for  which  we 
were  able  to  obtain  the  greatest  amount  of  data. 

Table  III.1  lists  some  of  the  physical  properties  of  samples  of 
type  T  and  samples  of  type  S.  Table  III. 2  lists  physical  properties  of 
samples  of  type  M.  The  dimensions  were  measured  using  scanning  electron 
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Figure  III. 11 .Schematic  diagram  of  the  experimental  set-up  used  to  measure  I-V  curves 


Figure  III.12.SEM  picture  of  burnt-out  microbridges.  The  fiber  in  (b)  is 
a  whisker  crystal  that  grew  out  of  the  film  which  had  been 
stored  for  more  than  a  year  before  this  picture  was  taken. 
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(III. 4) 
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CHAPTER  IV 


RESULTS  AND  DISCUSSION 

IV. 1. SAMPLES  OF  TYPE  S:  INTERACTING  PHASE-SLIP  CENTERS 

We  exanined  the  I-V  characteristics  of  each  of  the  five  segments  of 
the  bridge.  The  current  was  sequentially  applied  between  consecutive 
probes  and  the  voltage  measured  between  the  same  pair  of  probes  (for 
example,  C  to  L).  Figure  IV. 1  shows  a  set  of  I-V  curves  taken  in  this 
way  for  each  of  the  five  segments  of  sample  S-1 .  The  first  voltage  step 
in  each  of  the  I-V  curves  corresponds  to  the  appearance  of  a  phase-slip 
center  in  the  microbridge  and  the  other  voltage  steps  correspond  to  PSCs 
in  the  probes.  We  took  the  sane  set  of  I-V  curves  while  measuring  the 
voltage  between  the  ends  of  the  bridge  (leads  A  and  B)  in  order  to  check 
that  the  first  voltage  jumps  occur  in  the  bridge.  We  notice  that  the 
critical  current  of  the  probes  is  more  than  twice  the  critical  current 
of  the  microbridge,  as  expected  since  the  probes  are  about  twice  as  wide 
as  the  microbridge. 

For  each  segment  we  have  measured  the  critical  current  at  different 
temperatures  and  we  have  found  that  it  follows  the  relation 
IcoC(1-T/Tc)^2  predicted  by  the  Ginzburg-Landau  theory.  The 

2/g 

extrapolation  to  Ic=0  in  the  graph  of  I  versus  T  is  the  mean-field 
critical  temperature  T  .  The  results  (with  precision  of  +1mK)  for 
successive  segments  AC,  CL,  LD,  DM,  and  MB  of  the  sample  S-1  are  Tc= 
3.757,  3-745,  3.744,  3.746  and  3-746  K,  respectively.  In  the  same 
order,  the  slopes  of  Jc^3  versus  T  for  each  segment  are  1.9,  1.9,  2.2, 
2.3,  and  2.1  (yAA/(/<m)2)2/3(tnK)-1 . 

We  have  studied  the  interaction  of  a  PSC  in  one  segment  of  the 
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Figure  IV. 1.I-V  characteristics  of  each  one  of  the  five  segments  of 
sample  S-1 .  Variations  in  I  along  the  bridge  are  due  to 
imperfections  such  as  nonuniform  width.  The  inset  shows  the 
sample  geometry. 
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microbridge  with  PSCs  in  the  other  segments  of  the  same  microbridge. 
Consider,  for  example,  the  effects  that  a  PSC  in  segment  LD  has  on  the 
appearance  of  a  PSC  in  segment  CL.  We  measured  the  I-V  characteristic 
of  segment  CL  with  no  current  flowing  in  segment  LD  (figure  IV. 2,  curve 
I(LD)=0).  We  also  measured  the  same  I-V  for  fixed  values  of  I(LD)/0. 

No  effect  is  seen  as  long  as  I(LD)  is  less  than  the  critical  current  of 
the  segment  LD.  For  larger  values  of  I(LD)  we  see  a  change  in  the 
(apparent)  critical  current  of  CL.  The  critical  current  of  CL  is 
increased  or  decreased  depending  on  whether  the  current  in  the  two 
segments  is  flowing  in  the  same  direction  or  in  opposite  directions, 
respectively. 

We  notice  that  for  the  same  value  of  I(LD) ,  the  magnitude  of  the 
change  of  Ic(CL)  is  not  the  same  in  both  biasing  directions.  This  is 
because  the  heat  generated  in  LD  tends  to  decrease  the  critical  current 
of  CL  (by  increasing  its  temperature) ,  independent  of  the  relative 
direction  of  the  current.  (The  large  voltage  step  in  the  I-V  curves  for 
the  case  of  current  flowing  in  opposite  directions  result  from  the 
breakdown  of  superconductivity  in  probe  L,  where  I(LD)  and  I(CL)  add 
constructively) . 

In  the  sane  way  that  Jillie^  did  with  his  interacting  short 

microbridges,  we  decomposed  the  change  of  the  critical  current  of  each 

segment  A I  ,  into  a  symmetric  (even)  part,  A  I*,  and  an  antisymmetric 
c  c 

(odd)  part  Al^.  Here  the  symmetry  refers  to  the  effect  of  the 
relative  directions  of  the  current  in  the  source  and  detector  segments. 

We  define  In(0)  as  the  value  of  the  normal  current  at  the  center  of 
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opposite  directions,  respectively. 
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the  phase-slip  center  (source).  This  normal  current  is  a  fraction  of 
the  total  current  and  can  be  measured  directly  from  the  I-V  curves. 
Following  the  SBT  model^,  we  take  1^(0)  to  be  equal  to  the  voltage 
across  the  PSC  divided  by  the  differential  resistance  of  the  PSC  in  the 
linear  region.  In  segment  LD  we  have  found  that  In(0)=I(LD)-0.60Ic(LD) , 
which  is  very  close  to  the  SBT  prediction  In(0)=I-0.65Ic.  Similar 
results  were  found  in  other  segments. 

In  figure  IV. 3  we  have  plotted  AlS(CL)  and  Ala(CL)  as  a 

c  c 

function  of  In(0)  in  segment  LD.  Since  the  presence  of  a  PSC  in  CL 
reacts  back  on  the  PSC  in  LD,  for  each  point  in  figure  IV. 3  we  have 
determined  this  In(0)  from  the  voltage  measured  across  LD  just  at  the 
threshold  for  a  voltage  to  develop  in  the  I-V  curve  of  CL  at  the 
(shifted)  Ic(CL). 

First,  I  will  discuss  the  antisymmetric,  "heating- independent" 
part,  caused  by  the  quasiparticles  generated  in  the  phase-slip  center  in 
LD  that  diffuse  into  CL.  This  causes  the  magnitude  of  the  supercurrent 
to  change  in  order  to  maintain  the  total  current  through  CL  equal  to  the 
external  applied  current.  Consequently,  for  a  fixed  value  of  total 
current  through  CL  (before  the  appearance  of  a  phase-slip  center  in  CL), 
the  supercurrent  will  be  smaller  (or  larger)  than  the  current  applied, 
producing  an  increase  (or  decrease)  of  the  apparent  critical  current  of 
CL.  This  change  is  equal  in  magnitude  to  the  value  (measured  at  CL)  of 
the  quasiparticle  diffusive  current  generated  in  LD.  Therefore,  the 
value  of  the  slope  of  Ala(CL)  versus  1(0)  gives  the  fraction  f(x)  of 
the  quasiparticle  current  generated  in  LD  that  reaches  the  segment  CL. 
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Figure  IV. 3. Symmetric  and  antisymmetric  Al^  parts  of  the  change  of  the 

critical  current  of  CL  of  the  sanplecS-l  as  a  function  of  the  quasiparticle 
current  in  the  center  of  the  phase-slip  center  in  LD.  The  inset  shows  the 
linear  temperature  rise  at  CL  as  a  function  of  the  power  dissipated  in  LD. 
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In  a  uniform  filament,  this  fraction  should  be  equal  to  exp(-x/Ag#), 
where  x  is  the  distance  between  the  phase-slip  center  source  and  the 
detector.  Our  geometry  is  more  complicated  since  we  have  to  take  into 
account  the  effects  of  the  probes,  which  provide  another  channel  into 
which  the  quasiparticles  can  diffuse,  causing  a  more  rapid  decrease  in 
f(x)  frcm  one  observation  point  to  the  next.  The  effect  of  the 
phase-slip  center  in  DM  on  the  critical  current  of  CL  is  shown  in  figure 
IV. 4.  The  I-V  curves  of  CL  were  taken  for  different  fixed  values  of 
I(DM) . 

Figure  IV. 5  shows  the  results  for  f(x)  at  T=3-708fC.  By  definition, 
f(x=0)=1  at  the  sources,  and  it  decays  with  distance.  The  dotted 
straight  lines  are  guides  to  the  eye.  However,  their  parallelism  is 
consistent  with  the  assumptions  that  Aq*  is  about  the  same  everywhere 
along  the  bridge  and  that  the  same  fraction  of  quasiparticle  current  is 
lost  by  diffusion  into  each  of  the  four  probes  C,  L,  D,  and  M,  so  that 
f(x)  falls  in  geometric  progression  from  point  to  point.  The  horizontal 
error  bars  indicate  the  uncertainty  in  our  knowledge  of  the  exact 
location  of  the  core  of  the  phase-slip  center  in  each  segment. 

In  order  to  obtain  a  value  for  the  diffusion  length  Ag,  from  the 
data  shown  in  figure  IV. 5,  it  is  necessary  to  know  what  fraction  of 
quasiparticle  current  is  lost  to  the  probes.  This  can  be  determined 
theoretically  by  solving  the  steady-state  diffusion  equation  for  the 
charge  imbalance 


Q* 


(IV. 1) 
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This  diffusion  equation  for  Q*  is  the  dc  limit  of  the  more  general 
charge  imbalance  equation  derived  earlier1^. 

We  have  solved  this  diffusion  equation  numerically  for  the 
particular  geometry  of  sample  S-1  and  have  calculated  the  diffusion 
current  along  the  bridge,  normalized  to  its  value  at  the  center  of  the 
phase-slip  center.  We  solved  this  problem  in  the  case  when  the 
phase-slip  center  is  in  the  segment  LD  of  the  bridge  and  also  when  the 
phase-slip  center  is  in  the  segment  MB.  We  compared  the  results  for 
different  values  of  _Aq*  with  the  data  and  found  that  the  best  fit  is 
for  Aq*=(6.01  1.0)^im.  These  data  were  taken  at  T=3.708K, 
corresponding  to  1-T/Tc=0.01.  Using  the  literature  values®®, 

Vp=0. 65x10®  cm/sec  and®^  pH  =1 .05x10”®«tacm",  we  infer  that  for  sample 
S-1  the  mean  free  path  l=24x10~®  cm  and  that  the  measured  value  of 
at  T/Tc=0 . 99  corresponds  through  equation (II. 9) to 

2 

Xq,=  -^Q*  =  (71  2)x10~10sec 

VF^ 

Near  Tc,  '5q#  is  expected1 ®’ iy ’ ^  to  vary  as  ( 1-T/T^-1^2  ,  this  implies 
that  our  measured  value  corresponds  to 

Sq,  =  (0.7  +  0.2)x(1-T/Tc)_1/2x10'10  sec 
or 

Sq,  =  (1.2!  0.4)  A(0)/A(T)x10-10  sec 
Although  our  method  is  entirely  different,  this  result  is  in  excellent 
agreement  with  the  value  obtained  by  Clarke  and  Paterson®®  by  tunnel 
injection  technique. 
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We  turn  now  to  the  symmetric  change  of  the  critical  current  of  the 

detector  PSC,  to  confirm  that  it  can  be  accounted  for  by  the  heat 

generated  in  the  source  PSC.  We  can  estimate  the  temperature  rise  (AT) 

in  the  detector  for  a  given  power  generated  at  the  source  by  measuring 

A  Is  and  then  using  the  graph  I  (of  the  detector)  versus  T  to 
c  c 

determine  the  temperature  rise  in  the  segment  CL,  AT(CL),  caused  by  the 
phase-slip  center  in  the  segment  LD.  The  inset  in  figure  IV. 3  shows  the 
dependence  of  AT(CL)  on  the  power  dissipated  in  LD,  P(LD) .  The  data 
points  fall  on  a  straight  line,  the  slope  of  which  corresponds  to  a 
temperature  rise  in  CL  of  0.46mK  per  each  nW  of  power  dissipated  in  LD. 
The  curve  A 1^  versus  In(0)  in  figure  IV. 3  ,  calculated  using  this 
ratio,  gives  an  excellent  fit.  We  now  check  whether  this  fitted  value 
for  AT/P  is  reasonable. 

■? 

According  to  the  SBT  model  ,  the  temperature  rise  at  the  center  of 
the  PSC,  for  power  dissipation  IV,  is 


AT  .  1 _ _ 

IV  <2T\+AQ#)otw 


(IV. 2) 


i  /? 

where  r\  =  (Kd/<=0  is  the  thermal  healing  length,  K  is  the  thermal 

conductivity  of  the  metal,  o(  is  the  heat  transfer  coefficient  (to  the 
helium  bath  and  to  the  substrate)  per  unit  area,  and  w  and  d  are  the 
width  and  the  thickness  of  the  microbridge.  The  express  ion  (IV.  2)  was 
derived  for  a  one-dimensional  filament  and  can  not  be  compared  directly 
with  our  experimental  result  since  the  sample  S-1  has  a  more  complicated 
geometry . 
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In  order  to  estimate  the  magnitude  of  and  o<  ,  we  have  to  solve 
the  equation  governing  AT  in  the  sample  S-1  and  then  compare  with  our 
experimental  result.  Outside  the  source,  this  equation  has  the  same 
form  (see  reference  64)  as  the  diffusion  equation  for  Q*,  the  only 
change  being  replacement  of  Aq#  by  •  Therefore,  we  can  use  our 
numerical  solution  of  the  diffusion  equation  to  estimate  the  temperature 
rise  at  different  points  in  the  microbridge  for  different  values  of 
with  the  boundary  condition  that  dT/dx  =(2Kdw)“V(LD)  at  the  center  of 
the  PSC.  The  thermal  conductivity  can  be  calculated  from  p  by  use  of 
the  Wiederaann-Franz  law  K=  .  For  the  sample  S-1  , 

K=0 .22  watt/K  cm  at  TQ. 

We  found  that  our  temperature  rise  data  correspond  to  a  thermal 

healing  length  f[=(12i1)  ji m  and  a  heat  transfer  coefficient 

2 

0(=( 1.5 to. 2) watt/K  cm  .  This  value  is  consistent  with  the  experimental 

64  66 

values  obtained  by  Skocpol  ’  .  He  found  that  for  tin  films  on  glass 

or  sapphire  substrates  and  in  vacuun  or  in  direct  contact  with  the 

normal  helium  bath,  th-j  values  of  fall  in  the  range  1.0  to  3-5 
2 

watt/K  cm  . 

The  small  temperature  rise  in  the  microbridge  slightly  affects  the 
magnitude  of  the  antisymmetric  interaction,  since  Aq,  is  a 
temperature-dependent  parameter.  We  observed  this  effect  in  the  data 
points  of  Al*(CL)  in  figure  IV. 3.  Those  points  do  not  fall  exactly  on 
a  straight  line,  but  instead,  they  follow  a  curve  with  slightly 
increasing  slope  as  the  current  through  LD  increases.  This  variation  is 
small  and  it  is  reflected  in  the  error  bars  in  figure  IV. 5,  and  it  is 
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also  included  in  the  quoted  uncertainty  in  the  value  for  Aq*.  The 
straight  line  in  figure  IV. 3  is  an  average  of  all  points  including  the 
origin. 

Another  effect  that  should  be  taken  into  account  is  the  depression 

of  the  gap  due  to  the  presence  of  a  phase-slip  center,  as  observed 

41 

experimentally  by  Skocpol  o;id  Jackel  .  However,  this  effect  relaxes 
exponentially  within  a  coherence  length,  which  for  sample  S-1  is  equal 
to  1.3^.m  at  1=3.708  K.  If  we  assuned  that  the  gap  went  to  zero  at  the 
core  of  the  Fhase-slip  center  in  LD,  the  value  of  the  gap  at  the  middle 
of  CL  would  be  A  (x=7.8yKm)=  n-e"7'8/1 ,3)*0.998A-  and  the 

decrease  in  the  critical  current  of  CL  due  to  this  effect  alone  would  be 
less  than  IS.  This  is,  of  course,  very  small  compared  with  the  measured 
Al^(CL).  Therefore  we  can  conclude  that  the  heat  generated  in  the 
phase-slip  center  source  is  mainly  responsible  for  the  symmetric  change 
A Ig  in  the  detector . 
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IV. 2. SAMPLES  OF  TYPE  T: 

IV. 2.1. MEASUREMENTS  OF  THE  NONEQUILIBRIUM  POTENTIAL 
In  our  experiments  with  samples  of  type  T,  the  current  is  applied 
along  the  microbridge  and  the  voltage  is  measured  between  each  normal 
probe  and  one  end  of  the  microbridge.  When  the  current  exceeds  the 
critical  current  of  the  weakest  point,  a  phase-slip  center  is  formed. 

In  samples  with  a  narrow  segment,  like  the  one  shown  in  figure  III. 3, 
the  first  PSC  nucleates  in  that  segment.  In  the  set  of  I-V  curves  shown 
in  figure  IV.  6  we  can  see  that  the  biggest  voltage  drop  occurs  between 
probes  L  and  D. 

On  either  side  of  the  core  of  the  PSC,  where  it  makes  a  "step" 
change,  the  potential  of  the  pairs  is  constant  (see  figure  II. 1). 
Therefore,  the  voltage  measured  is  the  potential  difference  between  the 
pairs  and  the  quasiparticles  This  difference  is  related  to 

the  current  density  JR  in  the  following  way:  ?n  o<C^7  ( /^rf  • 

Equations  (11.11),  (11.12)  and  (II .13)  show  that  in  a  uniformly 
wide  filament,  Jn  and  (y*n-^<s)  decay  exponentially  with  the  same 
characteristic  length  A In  order  to  observe  this  exponential  decay, 
our  data  must  be  corrected  for  the  fact  that  the  middle  segment  of  the 
microbridge  is  narrower  than  the  rest.  Suppose,  for  example,  that  the 
width  of  the  middle  segment  is  equal  to  1/N  times  the  width  of  the  rest 
of  the  bridge;  then  for  a  given  current,  the  voltage  drop  in  that 
segment  would  be  N  times  larger  than  the  voltage  drop  in  the  same 
segment  if  the  bridge  were  uniformly  wide.  We  made  that  correction  in 
analyzing  the  data  from  samples  TN-1  and  TN-2.  Figure  IV. 7  shows  the 
results  for  sample  TN-2  before  and  after  such  correction  was  made.  We 


; 


i  is  applied  between  the 
Ltage  is  measured  between 
and  M  and  also  between  A  and  B. 
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chose  x=0  at  the  point  where  the  voltage  drop  is  half  the  total  voltage 
drop  between  the  ends  of  the  microbridge.  The  horizontal  error  bars 
indicate  the  size  of  the  normal  probes.  By  making  x^=0  in  equations 
(11.11)  and  (11.12),  we  can  write 


e  2 

(Ar/V  =  -  L  ex/Aj* 

e  2 


(x>0)  (IV. 3) 


(x<0)  (IV. 4) 


where  VQ  is  the  total  voltage  drop  between  the  ends  of  the  bridge.  To 
test  these  predictions,  in  figures  IV. 8  and  IV. 9  we  plot  log  v(x),  where 
v(x)=2|/<n-/Ksj /eVQ,  as  a  function  of  x,  for  samples  TN-2  and  TN-1, 
respectively.  The  straight  lines  show  that  the  decay  is  in  fact 
exponential.  The  slope  of  those  lines  is  a  direct  measure  of  the 
diffusion  length  of  the  quasiparticle  charge  Aq„ . 

In  sample  TU-1 ,  without  constriction,  the  first  phase-slip  center 
appeared  in  segment  MB.  Figure  IV. 10  shows  a  plot  of  log  v(x)  as  a 
function  of  x  for  x<0.  No  geometrical  correction  was  necessary  since 
TU-1  has  uniform  width. 

We  have  studied  the  temperature  dependence  of  Aq*  for  samples  of 
type  T  and  found  that  Aq#  diverges  at  the  critical  temperature.  We 
have  fitted  Aq#  to  the  expected  form  Aq,(T)=Aq,(0)  (1-T/Tc)_1/4  and 
the  corresponding  'Cq,  to  the  expression  'Sq#(T)  =  ’Sq*(0)(1-T/Tc)'"1/2. 

In  table  IV. I  a  list  is  shown  of  our  results  for  Aq*(0)  and  T>q#(0)  for 
our  samples.  In  order  to  facilitate  comparison  with  Clarke's  results, 
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we  have  also  quoted  our  results  in  the  form  A(t)/A(o,  using  the 
relation65  A(T)/A(0)X  1 .74(1  -T/Tc) 1/2 ,  near  Tc. 

The  values  of  ^(T)  A(T)/A (0)=1 .74  SQ,(0)  listed  in  table  IV. 1 
agree  extremely  well  with  the  value  obtained  for  the  sample  S-1 , 

A(T)/A  (0)=(1 .2  ± 0.4)x10-10  sec  and  with  the  value 
'Cq#A(T)/&  (0)  =  (1 .01  0.2)x10~10  sec  obtained  by  Clarke63  for  tin 
samples  very  close  to  Tc  in  experiments  where  a  nonequilibrium 
quasiparticle  population  had  been  created  by  current  injected  through  a 
tunnel  barrier. 

Very  close  to  T,,  can  be  written  as13,1^’15: 


->Q* 


:-!-BTqE  =  0.73 

TT  ACT) 


A  Co)' 

Act) 


(IV. 5) 


where  ^  is  the  inelastic  phonon  scattering  time  at  TQ.  Using  our 
experimental  results  from  all  our  samples,  our  value  for  G..  in  tin  is 
(1 .6- 0.4)x1Q~10  sec,  in  excellent  agreement  with  Clarke’s63  result, 
'5E=(1.i|i0.2)x10'10  sec. 

We  have  found  that  a  slightly  better  fit  of  our  data  is  obtained 
for  values  of  n  slightly  larger  than  1/4,  where  n  is  the  exponent  in  the 
expression  A^t(T)=  A'q»(0)  (1-T/Tc)"n.  The  results  for  n  and  jVq*(0) 
for  the  three  samples  are  listed  in  table  IV. 2.  In  figure  IV. 11  we  show 
the  values  of  Aq^T)  obtained  from  our  measurements  in  sample  TN-2.  We 
also  show  two  fits:  Ag,(T)=1 .3^.m(1-T/Tc)-0'26  plotted  as  a  full  line 
and  Aq,(T)=1 .5 /<m(1-T/Tc)"0,25  plotted  as  a  dashed  line.  Evidently, 
the  difference  between  the  two  fits  is  very  small.  For  samples  TN-1  and 
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TABLE  IV. 1. EXPERIMENTAL  VALUES  OF  QUASIPARTICLE  CHARGE  DIFFUSION  LENGTH 
AQ,  AND  RELAXATION  TIME  '&Q#. 


SAMPLE 

Aq#(0)^!11] 

'^Q*(0)[10"10sec] 

3Q.(T)A(T)/A(0)[10‘10sec] 

TN-1 

1  -3-  0.3 

0.710.3  ‘ 

1.21  0.5 

TN-2 

1.510.1 

0.61  0.1 

1.010.2 

TU-1 

1.210.2 

0.810.3 

1.410.6 

, 

TABLE  IV. 2. BEST  VALUES  OF  n  AND  A‘q»(0)  FOR 
\.(T)  =  A^(0)(1-T/Tc)-n. 

FIT  OF  OUR  DATA  TO 

SAMPLE 

n 

A’q*(0)[^  m] 

TN-1 

0.2710.04 

1.210.3 

TN-2 

0.2810.01 

1.310.1 

TU-1 

0.26  to. 04 

1.1  10.2 

—1  /u 

TU-1,  the  difference  between  the  best  fit  and  the  ( 1  — T/Tc ) 

divergence  is  even  smaller  than  for  the  sample  TN-2  and  the 

uncertainties  in  the  values  for  the  exponent  n  and  for  _A!q#(0)  are  also 

bigger  in  the  samples  TN-1  and  TU-1.  Therefore,  we  are  unable  to 

conclude,  from  these  results,  that  n  is  actually  larger  than  1/4. 

However,  if  the  difference  of  n  from  1/4  is  real,  it  might  be  a 

consequence  of  the  finite  size  and  dynamic  structure  of  the  core  of  the 

phase-slip  center  as  discussed  in  section  II. 4  and  also  of  the  elastic 

4-3 

processes  of  relaxation  J  of  charge  imbalance,  neither  of  which  is 
considered  in  the  derivation  of  equation  IV. 5.  According  to  Stuivinga 
et  al.  ,  in  the  presence  of  a  supercurrent  at  the  critical  value,  the 
relaxation  time  for  the  charge  imbalance  can  be  written  as 

'S0,  =  J1!!£c5-(1  +  1/(i(T))"1/2  (IV. 6) 

q  IT  ACT)  E  V 

where  (T)  is  defined  as 

(MT)  =  M  JL  _J - 1 -  (IV.7) 

'  8  Z-  k_T  ( 1-T/T  ) 

E  B  c  c 

Taking  this  effect  into  account,  the  diffusion  length  of  the 
quasiparticle  charge  will  diverge  at  Tc  according  to 

Aq#  oC  (1-T/Tc)“1/4(U1/^(T))"1/4  (IV. 8) 

From  our  experimental  results  for  Tc  and  ^g.,  we  found  n,  to  be 
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equal  to  1 .5x10”^(1-T/Tc)”^ .  In  figure  IV. 12  we  have  plotted 
[<1-T/TC)(1+1/(*>(T))]"1/4  and  (1-T/TC)"1/4  as  a  function  of  1-T/T,.  We 
can  see  that  the  stronger  temperature  dependence  of  the  theoretical  Aq# 
when  the  effects  of  the  supercurrent  are  included,  showing  how  this 
effect  can  be  partially  responsible  for  the  deviation  from  1/4  of  our 
experimental  values  of  the  exponent  n. 

If  we  use  equations  (IV. 6)  and  (IV. 7)  to  fit  our  data,  we  will  find 
that  the  best  value  for  is  slightly  larger  than  the  value  obtained 
using  equation  (IV. 5);  however,  this  difference  would  be  less  than  10% 
and  in  consequence,  we  take  'Sj.rO  .6  i  0.4)x10“10  sec,  obtained  from 
equation  (IV. 5),  to  be  the  correct  value  of  the  inelastic  scattering 
time  of  our  tin  samples.  It  is  also  important  to  point  out  that  in  the 
derivation  of  equation  (IV. 6)  it  is  assumed  that  the  supercurrent  is 
position- independent  which  is  not  true  in  a  phase-slip  center. 

IV. 2. 2. INJECTION  EXPERIMENTS 

We  have  injected  normal  current  into  the  microbridge  through  the 
normal  probes  in  samples  of  type  T.  We  measured  the  effect  of  that 
injected  current  on  the  appearance  of  the  first  PSC  in  the  microbridge. 
Figure  IV. 13  shows  a  set  of  I-V  curves  taken  in  sample  TN-1  for 
different  fixed  values  of  the  current  injected  through  D.  For 
I(DB)=10^.A  the  voltage  across  the  junction  D  is  260^. V;  therefore,  the 
energy  of  the  injected  quasiparticles  is  larger  than  the  energy  gap  in 
the  microbridge  at  that  temperature.  This  injected  current  has  the  same 
effect  on  the  critical  current  of  the  PSC  in  LD  as  the  effect  of  the 
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quasiparticle  current  generated  in  one  PSC  has  on  the  appearance  of 
other  PSCs  as  discussed  in  section  IV. 1. 

We  measured  the  change  in  the  critical  current  of  the  segment  LD 
when  the  normal  current  is  injected  through  D  and  also  when  the  normal 
current  is  injected  through  M.  We  extracted  the  part  that  depends  on 
the  relative  direction  of  the  currents,  Assuming  that  the 

nonequilibrium  current  decays  exponentially  between  M  and  D,  Aq#  is 
obtained  from: 

(°urrent  lnJected  thr°U8h  M)  -  vxpC-x^A,,)  (IV. 9) 
Ai  q(LD)  (current  injected  through  D) 

We  found  that  Aq#=(3.7£  1 m  at  1-T/TC=0.008  which  agrees  very 
well  with  the  value  Aq#=(4.3±  1 .0 )^im  that  we  obtained  from 
measurements  of  the  nonequilibrium  potential  created  by  a  PSC  in  sample 
TN-1 . 

In  order  to  write  equation  IV. 9  we  have  assumed  that  the  ratio  of 
quasiparticle  charge  injected  to  the  total  injected  current  is  the  same 
for  all  the  junctions.  This  assumption  is  valid  since  this  ratio 
approaches  unity  at  Tc,  independent  of  the  resistance  of  the  junction®®. 

We  have  also  measured  the  nonequilibrium  potential  induced  by  the 
injected  current.  The  current  is  applied  between  one  probe  and  one  end 
of  the  microbridge  and  the  voltage  measured  between  each  one  of  the 
other  probes  and  the  other  end  of  the  microbridge.  We  observed  the 
decay  of  this  potential  with  the  distance  from  the  injection  site  and 
assuming  that  this  decay  is  exponential,  the  magnitude  of  Aq»  is 
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obtained.  Within  the  limited  precision  of  this  measurements,  the 
results  in  samples  TN-1  and  TU-1  agree  with  the  values  of  Aq»  obtained 
from  measurements  of  the  nonequilibriun  potential  in  a  PSC.  These 
results  are: 


SAMPLE 

1-T/TC 

Ag#(from  injection  exp.) 

Ag*(from  PSC  data) 

TN-1 

0.0011 

(10*  4)^m 

(7.1  *  1.6)  yw.m 

TN-1 

0.0081 

(5*  2)^m 

(4.3i  1.0)  wm 

TU-1 

0.0071 

<5i  2)^m 

(4.0i0.7)y<m 
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IV. 3. SAMPLES  OF  TYPE  M:  EFFECTS  OF  MICROWAVES  ON  PHASE-SLIP 
CENTERS 

We  measured  the  I-V  curves  of  microbridges  of  various  lengths 
(samples  of  type  M)  at  various  temperatures  and  for  different  applied 
microwave  power.  Our  purpose  was  to  study  experimentally  the  effect  of 
the  length  of  the  microbridge  on  the  power  dependence  of  the  magnitude 
of  the  Josephson  steps  and  to  compare  with  the  Bessel-function  type  of 
dependence  for  short  microbridges. 

Josephson  steps  had  been  previously  observed  in  long  microbridges 

oft  llii  7 

of  tin  and  niobium  and  in  whiskers  of  tin  and  indium  .  However,  no 

quantitative  analysis  was  made  in  those  previous  works.  From  the 

theoretical  point  of  view,  we  are  still  lacking  of  a  model  that  predicts 

the  magnitude  of  the  current  steps  in  the  I-V  curves  of  superconducting 

filaments  irradiated  with  microwaves.  The  main  complication  has  to  do 

with  the  fact  that  the  oscillations  of  the  supercurrent  in  a  phase-slip 

center  are  not  restricted  to  a  point,  but  instead,  they  propagate  as 

damped  waves10  along  the  filament. 

44 

The  expression  (II. 36),  proposed  by  Tidecks  et  al.  ,  is  a 
simplification  of  the  problem  since  it  only  considers  the  time 
dependence  of  the  supercurrent.  Their  model  does  not  explain  the 
appearance  of  harmonic  and  subharmonic  current  steps. 

A  more  realistic  model  has  been  worked  out  by  Iv)ev  and 
Kopnin^®’^.  In  their  approximation,  discussed  in  section  II. 4,  the 
supercurrent  is  independent  of  x  for  x<x2(t),  where  x2(t)  oscillates 
between  0  and  XrDaxs:^r,^«  %  and  the  supercurrent  has  the  form  of  a 
saw-tooth  function  in  time.  For  (£5E>  ^>x>x2(t),the  oscillations  of 
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1  /o 

the  supercurrent  decay  with  x  until  they  are  negligible  for  E)  '  . 

r  is  the  pair  breaking  paraneter  and  it  is  taken  to  be  P <<1  in  the 
Ivlev  and  Kopnin  model,  i  g  is  the  penetration  length  of  the  electric 
field. 

The  magnitude  of  the  mfch  subharmonic  current  step  is  given  by 
equation  (11.30)  which  can  be  rewritten  in  the  form 

•  2  S4  V  eC  E~  (IV. 10) 

hxc  S" 

where  is  given  by  equation  (II. 31)  ,  is  the  normal  conductivity 
of  the  sample,  QT^  is  the  differential  conductivity  dl/dV,  and  E^,  is  the 
microwave  field  in  the  superconductor. 

According  to  equations  (IV. 10)  and  (11.36),  the  magnitude  of  the 
Josephson  steps  increases  indefinitely  with  the  microwave  power  applied. 
Experimentally,  heating  effects,  not  included  in  those  models,  reduce 
the  critical  current  of  the  microbridge  and  therefore,  reduce  the  width 
of  the  current  steps. 

In  the  experiments  that  I  am  reporting  here,  other  effects  were 
present  that  made  any  quantitative  analysis  of  the  power  dependence  of 
the  size  of  the  current  steps  very  difficult.  First  of  all,  the  samples 
of  type  M  exhibit  an  enhancement  of  the  critical  current  induced  by  the 

Cn  £Q 

microwave  radiation.  This  effect,  the  Wyatt-Dayem  effect0 '*  ,  had  been 

previously  observed  in  both  short  (L<  £(T))  and  long  (L»5  (T)) 
microbridges. 

Figures  IV. 14  and  IV. 15  show  sets  of  I -V  curves  of  samples  M-29 
ta’'en  at  diffi  -ent  powers  of  the  applied  microwave  radiation.  In  figure 


power.  See  text  for  discussion. 


Figure  IV.15.I-V  curves 
power.  See 
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IV.  14  we  can  see  the  first  Josephson  step  (V.j=21^iV)  at  35dB.  As  the 
power  is  increased,  the  critical  current  increases  and  the  voltage  jump 
becomes  larger  than  and  in  consequence,  the  first  Josephson  step  is 
not  seen.  Then,  as  the  microwave  power  is  increased  even  further,  the 
critical  current  is  depressed  and  the  current  steps  reappear.  The  noisy 
onset  of  the  voltage  in  the  I-V  curves  at  19dB  and  20dB  can  perbups  be 
explained  in  terms  of  the  variations  of  the  microwave  power  coupled  to 
the  sample  and  the  combination  of  two  opposite  effects:  the  enhancement 
of  the  critical  current  (Wyatt-Dayem  effect)  and  the  depression  of  the 
critical  current  caused  by  the  heat  generated  by  the  microwaves.  This 
noisy  onset  was  observed  in  the  narrowest  samples  where  heating  effects 
are  more  important. In  figure  IV. 15,  the  first  voltage  jump  is  larger 
than  V1  and  the  Josephson  steps  are  only  seen  for  levels  of  power  such 
that  the  critical  current  is  reduced  below  its  value  with  no  microwaves. 

The  I-V  curves  shown  in  figures  IV.  1*4  and  IV. 15  correspond  to  a 
sample  with  length  of  the  order  of  Aq«  and  therefore  only  one  PSC  can 
nucleate  there.  In  longer  microbridges,  where  two  or  more  PSCs  can 
coexist,  the  observation  of  microwave  steps  becomes  more  difficult.  On 
one  hand,  all  the  samples  exhibit  the  Wyatt-Dayem  effect  and  on  the 
other,  the  range  of  current  and  voltage  between  the  onset  of  the  first 
and  the  second  PSCs  is  usually  very  small. 

We  have  made  long  microbridges  with  a  narrow  segment  in  the  middle, 
using  the  mask  to  make  samples  of  type  T  but  doing  only  the  tin 
evaporation.  One  of  those  samples  in  shown  in  figure  III. 1(d).  In 
those  samples,  the  I-V  curves  show  a  long  range  of  constant  differential 
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resistance  after  the  first  voltage  jump.  We  observed  very  small  current 
steps  at  V1  in  one  of  those  samples. 

Our  data  are  summarized  in  table  III. 2  and  figure  IV. 16.  In  table 

III.  2,  the  order  of  the  observed  Josephson  steps  are  indicated.  Figure 

IV.  16  is  a  plot  of  the  length  versus  the  width  of  each  one  of  the 
samples  of  type  M,  normalized  to  the  temperature-dependent  G-L  coherence 
length  of  the  sample.  Each  sample  is  represented  by  a  line  instead  of  a 
point  because  the  microwave  experiments  were  performed  at  various 
temperatures  and  therefore,  at  various  values  of  P (T) .  The  long 
microbridges  with  a  narrow  segment  are  marked  with  horizontal  bars. 

The  thickest  lines  indicate  where  the  Josephson  steps  were 
observed.  Unfortunately,  we  have  not  been  able  to  identify  a  systematic 
explanation  for  the  observability  of  steps  in  only  a  fraction  of  samples 
of  similar  dimensions. 

The  vertical  dashed  line  indicates  the  critical  width  below  which 
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vortices  cannot  exist,  according  to  Likharev  ,  and  the  horizontal 
dashed  line  indicates  the  critical  length  above  which  the  current-phase 
relation  becomes  a  multivalued  function.  According  to  this 
classification,  the  upper  left  region  corresponds  to  long 
superconducting  filanents,  which  have  been  the  subject  of  this  thesis. 

No  conclusion  was  obtained  from  the  quantitative  analysis  of  the 
power  dependence  of  the  magnitude  of  the  steps  in  our  samples  because  of 
the  difficulties  discussed  above. 


Figure  IV. 16. Length  versus  width  of  sanples  of  type  M  normalized  to 
their  respective  temperature-dependent  coherence  length. 
The  thick  black  lines  indicate  where  Josephson  steps  were 
observed.  According  to  Likharev,  the  upper  left  region 
corresponds  to  one-dimensional  filaments. 


CHAPTER  V 


CONCLUDING  SUMMARY 

The  most  important  accomplishment  during  the  course  of  this  work 
has  been  the  experimental  determination  of  the  relaxation  time  of  the 
quasiparticle  charge  in  tin  from  measurements  of  the  one-diraensional 
diffusion  of  a  non-equilibrium  population  of  quasiparticles  created  by  a 
phase-slip  center  or  by  a  normal  current  externally  injected  in  a  long 
microbridge. 

In  one  experiment  we  have  observed  how  the  diffusive  currents  of 

quasiparticle  charge  and  energy  (heat)  generated  in  a  phase-slip  center 

(PSC)  in  a  tin  microbridge  affect  the  properties  of  other  PSCs  in  the 

same  microbridge.  We  have  measured  the  decay  with  distance  between  the 

two  PSCs  (the  source  and  the  detector)  of  the  quasiparticle  current,  and 

we  have  obtained  a  value  for  the  diffusion  length  of  Aq»=  (6  £  1  )^m  at 

T=0.99TC.  From  the  heating  effect  data  we  determined  that  the  heat 

2 

transfer  coefficient  is  o(  =  (1.5£0.2)watt/Kcnr  for  our  tin  samples  on 
glass  substrate  and  in  direct  contact  with  the  liquid  He  bath. 

In  another  experiment  we  have  measured  the  quasiparticle  potential 
in  tin  microbridges  with  attached  normal  tunnel  junction  probes.  We 
have  observed  the  exponential  decay  of  this  potential  with  distance  from 
the  center  of  the  PSC,  and  we  have  observed  the  divergence  at  Tc  of  the 
characteristic  decay  length.  We  found  that  j\gt(T)= Aq#(0) (1-T/Tc)”n 
with  ns 0.25. 

In  the  samples  with  attached  tunnel  junctions,  we  have  also 
measured  the  nonequilibrium  potential  induced  by  a  normal  current 
injected  through  one  of  the  junctions  and  we  observed  the  decay  of  this 
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potential  with  the  distance  from  the  injection  site.  The  effect  of  that 
injected  current  on  the  appearance  of  the  first  PSC  in  the  microbridge 
was  also  measured  and  the  decay  with  distance  between  the  injection  site 
and  the  PSC  was  observed.  From  those  injection  experiments  we  also 
obtained  a  value  of  Aq#  which  agrees  with  our  previous  results  on  the 
same  samples. 

The  relaxation  time  associated  with  the  experimental  value  of 
the  diffusion  length  Aq#  was  obtained  from  the  relation 

p 

-A.Q*/Vp l .  We  found  that,  for  all  our  samples  and  in  all  the 
experiments,  Sq#(T)=(0.7 t  0.2)x10“10sec(1-T/Tc)“1/2.  This  result 
indicates  that  the  relaxation  of  the  quasiparticle  charge  is  dominated 
by  electron-phonon  scattering  processes. 

The  inelastic  scattering  time  'Sj.  for  electrons  at  the  Fermi 

surface  and  at  T  was  determined  from  the  theoretical  relation 
c 

7SQ*(T)=4kgTc'Sj./nA(T) ,  which  can  be  written  as 

,5Q#(T)=0.42'SE(1-T/TC)"1/2.  We  found  that  ^=(1 .61  0.4)x10~10sec  in 
tin.  This  value  agrees  very  well  with  values  found  in  the  literature. 

Finally,  we  fabricated  tin  microbridges  of  various  lengths,  from 
L«:  1pm  to  L~  40 yu.m.  We  measured  the  I-V  curves  of  those  microbridges 
when  irradiated  with  microwaves,  and  we  observed  well-defined  Josephson 
steps  in  microbridges  up  to  10 ^m  long  and  tiny  steps  in  a  sample  36y^m 
long.  These  observations  prove  the  oscillatory  character  of  the 
supercurrent  at  the  core  of  the  PSC.  On  the  other  hand,  no  conclusion 
was  obtained  from  the  quantitative  analysis  of  the  power  dependence  of 
the  magnitude  of  the  steps  ,  because  of  the  presence  of  other  effects 


104 


such  as  the  Wyatt-Dayem  effect  and  heating  effects.  The  complications 
due  to  these  effects  also  prevented  us  from  being  able  to  predict  the 
observability  of  Josephson  steps  based  on  the  measured  physical 
properties  of  the  samples. 

In  summary,  we  have  performed  various  types  of  measurements  on 
phase-slip  centers  in  tin  microbridges  where  we  demonstrate  the 
exponential  decay  of  both  the  quasiparticle  current  and  the 
nonequilibriun  potential  and  the  existence  of  oscillations  of  the 
supercurrent  at  the  Josephson  frequency.  We  have  obtained  values  of  the 
relaxation  time  V  ,  the  magnitude  and  temperature  dependence  of  which 
agree  well  with  the  expected  relaxation  time  of  the  branch  imbalance 
mode  of  disequilibrium. 
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